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ANALYSIS  OF  COASTAL  SEDIMENT  TRANSPORT  PROCESSES  FROM 
WRIGHTSVILLE  BEACH  TO  FORT  FISHER,  NORTH  CAROLINA 


by 

T.C.  Winton,  I.B.  Chou,  G.M.  Powell,  and  J.D.  Crane 

I.  INTRODUCTION 

This  report  presents  a  comprehensive  engineering  analysis  of  the 
coastal  sediment  transport  processes  along  a  42-kilometer  segment  of  the 
North  Carolina  shoreline  from  Wrightsville  Beach  to  Fort  Fisher. 

Included  in  the  analysis  is  an  interpretation  of  all  available  data 
describing  the  littoral  processes,  longshore  transport,  and  the  behavior 
and  success  of  beach  nourishment  projects  at  Wrightsville  Beach  and  at 
Carolina  Beach,  North  Carolina. 

Several  coastal  engineering  studies  have  been  conducted  within  the 
study  area  to  assess  the  nearshore  coastal  processes  and  shoreline 
erosion  trends.  Vallianos  (1970)  investigated  the  influence  of  the 
manmade  Carolina  Beach  Inlet  on  the  volumetric  erosion  trends  of  the 
Masonboro  and  Carolina  beach  shorelines.  He  presented  a  preliminary 
assessment  of  the  impact  of  Masonboro  Inlet  north  jetty  on  the  longshore 
transport  trends  for  Wrightsville  and  Masonboro  beach  shorelines,  and  an 
evaluation  on  the  performance  of  the  1965  Carolina  Beach  beach  fill. 

Jarrett  (1977)  conducted  a  study  for  the  30-kilometer  segment  of 
shoreline  from  Wrightsville  Beach  to  Kure  Beach  in  relation  to  an 
environmental  assessment  of  coastal  erosion  as  affected  by  Carolina 
Beach  Inlet.  He  estimated  the  annual  rate  of  littoral  transport  between 
nine  littoral  cells  by  using  a  calibrated  energy  flux-wave  refraction 
sediment  budget  approach.  Jarrett  refined  Vallianos'  (1970)  bypassing 
rates  for  both  Masonboro  and  Carolina  Beach  Inlet  and  reassessed  the 
magnitude  of  the  impact  on  shore  process  of  manmade  changes  occurring 
during  the  study  period.  The  results  of  this  study  are  also  available 
in  reports  by  the  U.S.  Army  Engineer  District,  Wilmington  (1976;  1977). 

The  U.S.  Army  Engineer  District,  Wilmington  (1974),  presented 
historic  shoreline  changes  in  the  vicinity  of  Fort  Fisher  between  1865 
and  1973.  Several  plans  were  recommended  to  protect  the  historic  Fort 
Fisher  battlements  from  critical  dune  erosion. 

Large  quantities  of  data,  some  of  which  are  not  available  to  previous 
investigators,  were  evaluated  during  this  study.  Much  of  the  field  data 
were  collected  from  1964  to  1975  for  shoreline  erosion  studies  conducted 
by  the  U.S.  Army  Engineer  District,  Wilmington,  and  in  part  for  the 
Coastal  Engineering  Research  Center's  (CERC)  Beach  Evaluation  Program 
(BEP).  Profile  surveying  and  the  collection  of  other  data  used  in  this 
report  were  coordinated  by  CERC.  Data  evaluated  include  repetitive 
beach  profiles,  sand  data,  bathymetry  surveys,  wave  gage  records, 
dredging  records,  meteorological  records,  coastal  structure  design, 
coastal  georaorphological  studies,  shoreline  erosion  studies,  aerial 
photography,  and  beach  photography. 


Appendixes  A  to  G  present  a  graphic  description  of  the  shoreline 
changes  along  the  study  area  between  1964  and  1975.  These  plots  allow  a 
quantitative  assessment  and  interpretation  of  beach  response  to  seasonal 
climatic  changes,  storm  events,  beach-fill  projects,  and  coastal 
engineering  structures.  Long-term  trends  are  identified  and  used  to 
establish  a  sediment  budget  model  of  Wrightsville  and  Carolina  Beaches. 
The  analysis  of  the  excursion  distance  response  of  the  mean  low  water 
(MLW),  mean  sea  level  (MSL) ,  and  mean  high  water  (MHW)  contours  of 
profiles  along  Wrightsville  and  Carolina  Beaches  permitted  the  formu¬ 
lation  of  a  mathematical  description  of  post  beach-fill  performances. 

All  analyses  and  interpretations  of  results  are  included  in  this 
report.  Supplementary  data  are  provided  in  eight  unpublished  volumes 
(I  to  VIII)  which  are  available  from  the  CERC  technical  library. 

Volume  I  contains  five  sections:  Section  A  provides  a  beach  profile 
documentation  for  the  entire  study  shoreline;  Section  B  presents  storm 
histories  (accounts  of  the  major  storms  occurring  in  the  study  area); 
Section  C  provides  a  wave  refraction  analysis  of  the  area  including  wave 
gage  data  for  selected  wave  spectra  plots,  selected  data  from  CERC's 
Littoral  Environment  Observation  (LEO)  program,  and  wave  refraction 
plots;  Section  D  presents  plots  and  tabulated  values  of  the  gross 
northerly  and  southerly,  and  the  net  longshore  energy  flux  distribution; 
and  Section  E  provides  data  on  volumetric  changes  which  occurred  within 
all  inlets  along  the  study  area.  Comparative  short  and  long  beach 
profiles,  beach  profile  data,  MSL  excursion  rate  tables,  MSL  volumetric 
change  plots  and  tables,  and  selected  sand  data  are  presented  for 
Wrightsville  Beach  (Vols.  II,  III,  and  IV),  Masonboro  Beach  (Vol.  V), 
Carolina  Beach  (Vols.  VI  and  VII),  Kure  Beach  (Vol.  VIII,  Sec.  I),  and 
Fort  Fisher  (Vol.  VIII,  Sec.  J). 


II.  STUDY  AREA 

The  study  area  is  part  of  the  tidewater  region  of  the  Atlantic 
Coastal  Plain,  consisting  of  a  series  of  low,  narrow,  sandy  barrier 
islands  and  peninsular  beaches  located  in  New  Hanover  County,  North 
Carolina.  The  islands  front  the  Atlantic  Ocean  just  north  of  Cape  Fear 
and  are  separated  from  the  mainland  by  either  the  Cape  Fear  River 
estuary  or  by  Myrtle  Grove,  Masonboro,  Greenville,  and  Middle  Sounds. 
The  five  coastal  sites  in  the  42-kilometer  study  are  (from  north  to 
south)  Wrightsville  Beach,  Masonboro  Beach,  Carolina  Beach,  Kure  Beach, 
and  Fort  Fisher.  Figure  1  shows  the  study  area  and  the  location  of  the 
five  study  segments. 


The  beach  sands  are  generally  fine  and  composed  of  quartz  sand  with 
a  shell  content  ranging  from  0  to  42  percent.  The  direct  sources  of 
littoral  materials  for  the  study  area  are  the  adjacent  beaches,  dunes, 
and  bluffs  (direction  of  transport  depending  on  direction  of  wave 
attack)  as  a  result  of  erosion,  and  the  nearshore  ocean  bottom  areas, 
from  which  material  is  brought  onto  shore.  A  complete  description  of 
the  geomorphology  and  geologic  history  of  the  study  area  has  been 
summarized  by  Pierce  (1970). 
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Based  on  data  recorded  by  CERC's  wave  gage  located  at  Wrightsville 
Beach,  the  annual  significant  wave  height  is  0.76  meter  (2.5  feet). 

Wave  observations  along  Wrightsville  Beach  indicate  that  98  percent  of 
the  observed  wave  energy  approaches  from  the  eastern  and  southeastern 
quadrants.  The  dominant  direction  of  littoral  transport  is  from  north 
to  south;  however,  reversals  in  transport  direction  along  the  beaches  do 
occur.  The  mean  and  spring  tidal  ranges  are  1.2  and  1.4  meters, 
respectively;  the  difference  between  MSL  and  MLW  is  0.57  meter. 

Wrightsville  Beach  is  about  6.75  kilometers  in  length,  with  an 
average  dune  height  of  4  meters  above  MSL.  The  beach  faces  approxi¬ 
mately  east-southeast,  has  an  average  beach  slope  from  MHW  to  the 
-6.0  meter  (MSL)  depth  contour  of  1  on  37.2,  and  contains  beach  sedi¬ 
ments  with  a  mean  grain  size  of  0.27  millimeter.  The  ocean  shoreline  of 
Wrightsville  Beach  was  modified  in  1965  by  the  construction  of  a  hurri¬ 
cane  and  storm  protection  project.  Initially,  2,288,000  cubic  meters  of 
fill  material  was  placed  along  5,100  meters  of  beach  north  of  Masonboro 
Inlet  with  artificial  dune  heights  constructed  to  an  approximate  eleva¬ 
tion  of  +2.5  meters  (MSL)  for  storm  protection  purposes.  The  northern 
transition  section  included  the  closure  of  Moore  Inlet,  which  had 
previously  separated  Wrightsville  Beach  from  Shell  Island.  In  spring 

1966,  an  additional  244,000  cubic  meters  of  fill  material  from  the 
Masonboro  Inlet  was  placed  between  Johnnie  Mercer's  Pier  and  Crystal 
Pier.  In  October  1966,  a  final  deposition  of  32,100  cubic  meters  of 
material  from  the  estuarial  area  behind  Shell  Island  was  placed  along 
the  northernmost  610  meters  within  the  town  limits  of  the  Wrightsville 
Beach  project  shoreline. 

In  1970,  a  renour ishraent  of  the  central  shoreline  of  Wrightsville 
Beach  was  required.  A  total  of  1,053,600  cubic  meters  of  fill  material 
obtained  from  a  shoal  in  the  Banks  Channel  and  the  sound  area  behind 
Shell  Island  was  placed  on  the  beach,  beginning  at  a  point  approximately 
1.83  kilometers  north  of  Masonboro  Inlet  and  extending  to  the  northern 
city  limits  of  Wrightsville  Beach.  Figure  2  is  an  aerial  photo  strip 
map  showing  the  Wrightsville  Beach  shoreline. 

Masonboro  Island  is  bordered  by  Masonboro  Inlet  to  the  north,  and  by 
Carolina  Beach  Inlet  (opened  in  1952  by  local  interest  groups)  to  the 
south  (Figs.  2  and  3).  It  is  a  very  narrow,  low-lying  uninhabited 
island  approximately  12.5  kilometers  long  with  a  shoreline  orientation 
from  north-northeast  to  south-southwest.  The  natural  dune  heights  along 
the  island  range  from  3  to  10  meters  (MSL),  and  the  median  grain  size  is 
0.34  millimeter.  The  average  beach  slope  is  approximately  1  on  59. 

Carolina  Beach  is  located  just  south  of  the  Carolina  Beach  Inlet  and 
extends  about  4.3  kilometers  southward  to  Kure  Beach  (Figs.  3  and  4). 

The  northern  end  of  Carolina  Beach  has  experienced  high  erosion  rates 
since  the  opening  of  Carolina  Beach  Inlet  (Vallianos,  1970),  which  have 
affected  the  efficiency  of  a  hurricane  and  shore  protection  project 
constructed  in  1965.  The  4.27  kilometers  of  shoreline  fronting  the  town 
of  Carolina  Beach  was  nourished  with  about  2,014,000  cubic  meters  of 
fill  material  obtained  from  the  Carolina  Beach  harbor.  However,  by 

1967,  erosion  of  the  northern  1.2  kilometers  of  the  project  beach  was  so 
severe  that  emergency  action  was  required.  Approximately  314,000  cubic 
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Figure  2.  Aerial  photo  map  of  study  area  from  Figure  Eight  Island 
to  Masonboro  Beach,  North  Carolina. 
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Figure  3.  Aerial  photo  map  of  study  area  from  Masonboro  Beach 
to  Carolina  Beach,  North  Carolina. 
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Figure  4.  Aerial  photo  map  of  study  area  from  Carolina  Beach  to  Fort  Fisher, 
North  Carolina. 
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meters  of  fill  material  was  distributed  there  and  83,400  cubic  meters  of 
sand  wa s  placed  to  form  a  new  520-meter  transition  section  from  the 
original  northern  limits  of  the  project  beach.  A  temporary  wooden  groin 
was  constructed  at  the  transition  junction  between  the  two  fill  sites. 

Despite  the  1967  emergency  action,  serious  erosion  continued, 
requiring  the  supplemental  emergency  construction  in  1970  of  a  335-meter 
rubble-mound  seawall  extending  southward  from  the  northern  boundary  of 
the  project.  In  conjunction  with  the  seawall  construction,  264,500 
cubic  meters  of  fill  material  from  the  sediment  trap  located  inside 
Carolina  Beach  Inlet  was  placed  along  the  northern  1.2  kilometers  of 
shoreline.  By  late  spring  1971,  the  southern  3.47  kilometers  of  the 
project  beach  had  been  partially  restored  with  approximately  581,000 
cubic  meters  of  material  from  a  borrow  area  located  in  the  Cape  Fear 
River.  The  rubble-mound  seawall  was  extended  an  additional  290  meters 
southward  in  1973.  The  severe  erosion  trend  of  the  northern  project 
limits  continued  despite  the  numerous  remedial  measures  taken. 

Kure  Beach  has  a  shoreline  about  4.25  kilometers  in  length,  and  is 
situated  between  Carolina  Beach  to  the  north  and  Fort  Fisher  to  the 
south  (Fig.  4).  The  city  of  Kure  Beach  and  the  unincorporated  towns  of 
Wilmington  Beach  and  Hanby  Beach  are  located  in  this  segment.  Dune 
heights  average  2.5  meters  above  MSL  along  this  segment;  beaches  have  a 
median  sand  grain  size  of  0.30  millimeter  and  an  average  beach  profile 
slope  of  1  on  30.  The  beaches  along  this  shoreline  remained  relatively 
stable  during  the  study  period. 

Fort  Fisher,  the  southernmost  segment  of  shoreline  studied,  is 
approximately  6.25  kilometers  long  and  extends  southward  from  Kure  Beach 
to  just  north  of  New  Inlet  (Figs.  4  and  5).  The  mean  grain  size  of  the 
beach  sand  is  0.27  millimeter  and  the  average  slope  is  approximately  1 
on  36.  The  northern  1.6  kilometers  of  shoreline  is  a  sandy  beach,  mostly 
undeveloped,  which  varies  in  width  from  27  to  55  meters.  This  section 
remained  relatively  stable  during  the  study  period.  The  central  stretch 
of  beach  contains  the  historic  remains  of  a  Confederate  Army 
fortification  known  as  Fort  Fisher,  which  was  built  adjacent  to  New 
Inlet.  Since  the  closure  of  this  inlet  in  1883,  rapid  erosion  exposed 
an  outcrop  of  coquina  rock  located  adjacent  to  the  remains  of  the  fort 
(Fig.  1).  The  sandy  beach  fronting  Fort  Fisher  varied  in  width  from  0 
to  45  meters  during  mean  tide  levels,  and  the  sand  bluff  along  the 
backshore  continued  to  erode  at  a  critical  rate,  thus  requiring  con¬ 
struction  of  an  emergency  rubble  revetment.  In  July  1965,  additional 
rubble  was  placed  along  both  the  northern  and  southern  flanks;  11,500 
cubic  meters  of  sand  was  also  placed  along  213  meters  of  shore  north  of 
the  revetment.  In  May  1967,  an  extratropical  cyclone  caused  severe 
erosion  to  the  1965  emergency  fill  which  required  placement  of  another 
11,500  cubic  meters  of  sand  along  the  same  beach  section.  In  1970, 
further  emergency  measures  were  implemented  by  placement  of  a  limestone 
revetment  along  a  part  of  the  upland  bluff  which  had  previously  been 
protected  by  the  beach  fills.  The  southernmost  4.58  kilometers  of  shore 
is  an  accreting  sandspit  characterized  by  low  topography  and  a  sandy 
beach  with  widths  between  60  and  275  meters. 

The  study  area  and  the  beach-fill  projects  are  further  described  in 
Vallianos  (1970),  U.S,  Army  Engineer  District,  Wilmington  (1970,  1974, 
and  1977),  and  Jarrett  (1977). 
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Figure  5.  Aerial  photo  map  of  study  area  from  Fori  Fisher 
to  Cape  Fear,  North  Carolina. 
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III.  DATA  COLLECTION 


1 .  Beach  Profiles. 

The  sediment  budget  analysis  performed  in  the  study  area  was  based 
on  the  beach  profile  data  provided  by  CERC.  Beach  surveys  were  taken 
at  241  stations  along  the  shoreline,  and  each  profile  was  perpendicular 
to  the  local  shoreline.  The  survey  stations  were  numbered  sequentially 
from  north  to  south  and  were  prefixed  by  the  abbreviation  of  the 
corresponding  beach  nane;  e.g.,  WB  for  Wrightsville  Beach  (50  stations), 
MB  for  Masonboro  Beach  (31  stations),  CB  for  Carolina  Beach 
(119  stations),  KB  for  Kure  Beach  (20  stations),  and  FB  for  Fort  Fisher 
Beach  (21  stations).  Station  CB2  would  therefore  represent  the  second 
station  from  the  north  in  Carolina  Beach.  Figures  6  and  7  show  the 
relative  locations  of  all  the  stations. 


The  beach  surveys  were  conducted  by  contractor  for  U.S.  Army 
Engineer  District,  Wilmington,  from  1963  to  1975.  Most  profiles  were 
measured  by  level  and  tape  and  extended  to  only  about  2.4  meters 
(8  feet)  or  less  below  MSL.  These  profiles  were  referred  to  as  short 
profiles.  Long  profiles  were  measured  to  a  depth  of  12.2  meters 
(40  feet)  using  a  depth  sounder.  Table  1  shows  the  survey  stations, 
along  with  CERC's  station  reference  codes,  which  indicates  long  profiles 
by  the  letter  L. 


About  2,952  repetitive  beach  profiles  were  taken  during  399  surveys, 
including  2,815  short  profiles  and  137  long  profiles.  Table  2  shows  the 
number  of  short  and  long  profiles  for  each  beach.  Table  2  and  Figures  6 
and  7  show  that  Wrightsville  and  Carolina  Beaches  have  much  better 
temporal  and  spatial  resolution  than  the  rest  of  the  study  area.  Of  the 
entire  beach  data,  89  percent  of  the  profiles  were  taken  on  Wrightsville 
and  Carolina  Beaches.  The  Fort  Fisher  Beach,  Kure  Beach,  and  Masonboro 
Beach  profile  data  were  of  insufficient  quantity  to  permit  a  valid 
analysis. 


All  data  are  available  in  supplementary  data  Volumes  I  to  VIII  from 
the  CERC  library. 


2.  Wave  Data. 


The  wave  climate  data  for  the  study  area  are  from  the  following 
sources: 

(a)  A  CERC  wave  gage,  located  on  Johnnie  Mercer's  Pier  at 

Wrightsville  Beach,  which  operated  from  March  1971  to  February 
1975.  The  gage  was  located  in  5.2  meters  (17  feet)  of  water, 
and  the  recorded  wave  data  represent  approximately  all  waves 
reaching  Wrightsville  Beach  from  all  seaward  directions. 
However,  wave  direction  could  not  be  differentiated  by  the 
gage.  The  wave  gage  data  for  this  study  with  selected  wave 
spectral  plots  are  presented  in  supplementary  data  Volume  I. 
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Figure  6.  Profile  station  location  map,  WB1  to  MB27. 
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Table  1.  Cross  references  for  beach  profile  data 


Note — Coastal  structures  at  profiles  WB17  (Johnnie  Mercer's  Pier) 
WB45  (Crystal  Pier),  CB42  (groin),  and  CB57  (fishing  pier). 
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Frying  Pan  Shoals  Light  Tower.  The  wave  data  with  the  monthly 
wave  statistics  were  provided  by  CERC. 


(c)  Long-term  deepwater  wave  statistics  provided  in  the  Summary  of 
Synoptic  Meteorological  Observations  (SSMO)  (U.S.  Naval  Weather 
Service  Command,  1975). 


CERC's  wave  observation  program  at  Wrightsville  Beach  provided 
visual  observations  of  wave  conditions,  recorded  daily  at 
Johnnie  Mercer's  Pier  between  June  1970  and  December  1973. 

CERC  provided  the  monthly  statistical  analysis  of  these 
shore-based  wave  observations  including  breaking  wave  height, 
period,  and  direction.  The  wave  data  collected  at  Wrightsville 
Beach  during  the  study  period  are  available  in  supplementary 
data  Volumes  II,  III,  and  IV. 


3 .  Beach  Sand  Data. 

Beach  sand  data  for  certain  profiles  within  the  study  area  from  1969 
to  1971  were  provided  by  CERC.  Samples  were  collected  along  the  profile 
azimuth  from  the  dune  crest,  the  berm,  and  at  MHW ,  MSL,  MLW,  -1.8  meters 
(-6  feet)  (MLW),  -3.66  meters  (-12  feet)  (MLW),  and  -5.49  meters 
(-18  feet)  (MLW).  Frequency  of  sand  sample  collection  was  not  con¬ 
sistent  from  beach  to  beach  or  from  profile  to  profile.  The  sand  was 
analyzed  for  basic  engineering  properties  including  grain-size  distribu¬ 
tion,  median  grain  size,  standard  deviation,  fall  velocity,  and  compo¬ 
sition.  Grain-size  analyses  are  summarized  in  Table  3.  The  complete 
sand  data  are  presented  in  supplementary  data  volumes  for  each  beach 
segment  (except  for  Kure  Beach). 
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NOTE — t p -  -lc^D,  where  D  =  sand  diameter  in  millimeters. 
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IV.  ANALYSIS  OF  BEACH  PROFILE  DATA 


1.  Excursion  Distance  Technique. 


If  successive  aerial  photos  of  a  beach  face  are  compared  with  each 
other  and  a  change  in  location  of  the  beach  is  noted,  then  this  change 
is  indicative  of  either  a  period  of  erosion  or  accretion.  Horizontal 
displacement  of  the  planform  position  of  any  one  point  on  the  beach, 
from  one  survey  to  another,  is  the  excursion  distance  for  that  point  for 
the  survey  period.  On  an  accreting  beach,  the  excursion  distance  of  a 
point  relative  to  its  initial  position  is  positive,  and  on  an  eroding 
beach,  it  is  negative.  The  rate  of  change  of  the  excursion  distance 
with  time  is  the  excursion  rate. 


If  successive  beach  profiles  are  reduced  to  a  common  base  line,  the 
excursion  distance  of  each  point  on  the  profile  indicates  the  magnitude 
of  the  onshore-offshore  movement.  The  relative  magnitude  of  the  excur¬ 
sion  distances  between  two  or  more  points  on  the  same  profile  identifies 
and  quantifies  the  change  in  beach  slope  between  those  points.  Beach 
excursions  can  be  converted  to  volumetric  changes  for  the  entire  active 
profile  by  applying  to  the  excursion  distances  a  volumetric  equivalent 
factor.  This  factor  was  developed  from  measured  changes  at  two  piers 
located  along  Wrightsville  Beach  (U.S.  Army  Engineer  District,  Wil¬ 
mington,  1977),  which  showed  that  for  a  closure  depth  of  approximately 
8.23  meters,  each  meter  of  excursion  was  equivalent  to  8.23  cubic  meters 
of  change  for  the  entire  active  profile  per  meter  of  beach  front. 
Equivalently  in  English  units,  for  a  closure  depth  of  27  feet,  each  foot 
of  excursion  was  equivalent  to  1  cubic  yard  of  change  for  the  entire 
active  profile  per  foot  of  beach  front.  Consequently,  excursion 
distance  analysis  is  a  simple  but  powerful  technique  which  is  used  to 
identify  and  quantify  both  long-term  beach  changes  and  the  response  of  a 
beach  to  short-term  impacts  resulting  from  storm  activity,  beach  fills, 
and  other  man- induced  changes. 


Meaningful  interpretation  of  excursion  distance  plots  can  only  be 
performed  if  known  short-term  or  sudden  impact  events  are  identified  and 
accounted  for  within  the  analysis.  In  order  to  do  this,  all  major 
erosion-causing  storms  and  all  man-related  activities  which  cause 
erosion-accretion  during  the  study  period  must  be  abstracted  from  the 
historical  records. 


Table  4  lists  all  beach-fill  changes  reported  along  the  study  area 
beaches  from  1965  to  1974.  The  initial  fill  excursion  distances  in  the 
table  were  estimated  by  applying  the  volumetric  equivalent  factor  of 
8.23  cubic  meters  of  change  for  each  meter  excursion  per  meter  beach 
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Beach-fill  evaluation 


front.  The  excursion  loss  due  to  sorting  was  determined  in  the  same 
manner  from  estimates  on  the  volume  of  beach  fill  lost  due  to  sorting  or 
from  volumetric  loss  calculations  based  on  the  critical  ratios  of  the 
beach-fill  material.  Note  that  the  initial  fill  distance,  the  excursion 
loss  (due  to  sorting),  and  the  net  fill  excursion  are  only  comparison 
estimates  and  should  not  be  considered  as  absolute  values.  Figure  8 
shows  the  spatial  distribution  of  the  beach-fill  excursions  along  the 
study  area,  with  an  obvious  concentration  of  fill  activity  in  front  of 
the  townships  of  Wrightsville  Beach  and  Carolina  Beach.  Areas  of 
reported  net  beach  fill  are  shown  to  extend  in  some  places  to 
approximately  100  meters  (300  feet).  Because  these  values  only  reflect 
the  fill  excursion  remaining  after  the  initial  loss  period  and  do  not 
consider  the  fill  loss  due  to  storm-induced  or  long-term  (annual) 
erosion  rates,  they  are  slightly  misleading.  Most  fills  were  placed 
after  the  previous  fill  had  been  severely  eroded  away. 

Table  5  presents  all  historical  events  influencing  beach  volumes 
since  1965,  with  a  brief  description  of  each  event.  Storms  were 
included  in  this  table  only  if  noted  beach  erosion  occurred,  if  asso¬ 
ciated  storm  surge  was  noted,  or  if  the  windspeeds  were  in  excess  of 
80  kilometers  per  hour  (50  miles  per  hour).  A  complete  list  of  all 
storms  during  the  study  period  is  available  in  supplemental  data 
Volume  I,  Section  B. 

3 .  Excursion  Distance  Analysis. 

Selected  beach  profiles  from  all  stations  were  plotted-  at  a  small 
scale  and  visually  checked  for  accuracy  and  acceptability  of  data 
points.  Larger  scale  profiles  were  then  drawn  to  compare  sequential 
outlines.  Areas  of  erosion  from  one  sequential  profile  to  the  next  were 
highlighted  by  a  dot-screen  pattern.  Typical  short  and  long  beach 
profile  plots  are  shown  in  Figures  9  and  10.  All  of  the  larger  scale 
plots  of  the  short  and  long  beach  profiles  are  contained  within  sup¬ 
plemental  data  Volumes  II  to  VIII. 

A  common  base  line  was  established  for  each  sequential  profile  and 
the  horizontal  distance  from  that  base  line  to  the  location  of  the  MHW, 
MSL,  MLW,  -1.83  meters  (-6  feet),  -3.66  meters  (-12  feet),  and 
-5.49  meters  (-18  feet)  contours  were  calculated.  These  distances  were 
plotted  against  time  of  measurement,  and  the  relative  distance  between 
the  first  and  subsequent  distances  represents  the  excursion  distance 
through  time  for  each  contour. 

A  sample  plot  from  each  beach  is  shown  in  Figures  11  to  15.  A 
linear  regression  ("least  squares")  line  which  mathematically  "best 
fits"  all  data  points  is  drawn  on  these  plots.  One  straight  line  is  not 
representative  of  the  average  excursion  rates  between  the  years  1965  and 
1975,  especially  for  Wrightsville  and  Carolina  Beaches. 

When  few  data  points  exist,  the  scatter  due  to  seasonal  fluctua¬ 
tions,  prior  storm  erosion,  etc.,  can  totally  mask  the  longer  term  or 
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Figure  8.  Distribution  of  (reach  Alls  along  study  area. 


Table  5.  Historical  events  affecting  beach  volumes  during 


study  period, 1965-1975. 


Data 

Locat ion 

Cantata 

rm 

Wright  evil  la  Beach 

Beach  fill  1. 9-7.0  ka;  47-*  nat  ex cur a ion 

Carolina  Baach 

Baach  fill  22.2-26.5  ka;  32-*  nat  axcuraion 

■I  ynn  w T 

Storm;  high  wind,  rain,  baach  eroaion 

MM: 

Fort  Fiahar  Baach 

Baach  fill  and  ravetaent  32.7-33.0  ka; 

6.5-n  nat  axcuraion 

1966 

Spring 

Wrightaville  Baach 

Beach  fill  3.4-6. 1  ka;  10-a  net  axcuraion 

Spring 

Maaonboro  Inlet 

Completion  of  Maaonboro  jetty 

10-11  Juna 

Tropical  Stora  Alaa  paaaad  offahore 

9  July 

Wrightaville  Baach 

Stona;  147-ka/h  (92  ai/h)  vinda 

Oct. 

Beach  fill  3. 4-4.0  ka;  6.5-a  axcuraion 

1967 

Mar. 

Carolina  Beach 

Beech  fill  21.7-23.5  ka;  10.5  a  n«t 
excursion 

19  Mar. 

Stora;  71-112  ka/h  (45  70  ai/h)  vinda 

29  May 

Bxtratropical  cyclone;  aavere  eroaion 

Oct. 

Fort  Fiahar  Baach 

Beach  fill  32.7-33.0  ka:  6.5-a  axcuraion 
Stora;  96-ka/h  (60  ai/h;  vinda 

24  Mov. 

28  Dec. 

Stora;  122-ka/h  (76  ai/h)  vinda 

1968 

7-12  Juna 

Tropical  Stora  Abby 

Aug. 

Carolina  Beach 

Beach  fill  23.0-23.7  ka;  13-a  net  axcuraion 

19-20  Oct. 

Hurricane  Cladya 

1969 

1-2  Nov. 

Stora;  96-ka/h  (60  ai/h)  vinda 

1970 

Mar. -May 

Wrightaville  Beach 

Bench  fill  2. 7-4. 6  ka;  epprox.  31.5-a  net 
excursion 

16-17  Aug. 

Stora;  2.5-a  (8  ft)  vavea,  riptide.; 

112-ka/h  (70  ai/h)  vinda 

30-31  Oct. 

Stora;  beach  eroaion 

Dec. 

Carolina  Beach 

Beach  fill  22.2-23.5  ka;  21-a  net  axcuraion 

Dec . 

Carolina  Beach 

Coaptation  of  rubble-aound  aeavall 

Fort  Fiahar  Beach 

Liaaatone  ravetaent  added 

1971 

26-30  Jan. 

Stora;  near  hurricane-force  vinda 

13  Feb. 

Stora;  near  hurricane-force  vinda 

5-7  Apr. 

Stora;  109-ka/h  (68  mi/h)  vinda 

Mar. 

Carolina  Beach 

Beech  fill  epprox.  23.0*26.5  ka;  11-a  net 
excursion 

16-18  Aug. 

Stora;  3-a  (10  ft)  aeaa 

27  Aug. 

Tropical  Stora  Dora;  96-ka/h  (60  ai/h) 
vinda,  1.2-a  (4  ft)  aurge 

Oct. 

Hurricane  Ginger;  147-ka/h  (92  ai/h)  vinda, 
1.2-a  (4  ft;  aurge 

1972 

24  July 

Stora;  83-kn/h  (52  ai/h)  vinda 

1973 

9-10  Fab. 

Stora;  80-ka/h  (50  ai/h)  vinda,  high  sees, 
erosion 

22  Mar. 

Stora;  3-4-a  (10-12  ft)  aeaa,  high  eroaion 

Sept . 

Carolina  Baach 

Ext ana ion  of  rubble-aound  aeavall 

1974 

30  Nov .-1  Di 

e. 

Stora;  eroaion 

MOTE:  Data*  of  beach  fill*,  coaatal  eonatruction,  ate.  ara  fivan  only  aa 
aonth  or  aaaaon  in  ich  they  were  coapleted.  Data*  of  itom  ara  given  aa 
calendar  data. 
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Figure  14.  Distance  from  the  base  Una  to  stated  contours  at  KB  17, 
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Figure  15.  Distance  from  the  base  line  to  stated  contours  at  FB  10, 


man-influenced  excursion  rates.  The  plots  (Figs.  11  to  14)  are,  in  one 
way,  atypical  of  all  profile  plots  taken  along  each  beach  because  each 
of  these  profiles  has  some  data  taken  below  MLW,  whereas  the  majority  of 
profiles  along  the  entire  study  area  do  not.  This  means  that  analysis 
of  contours  below  MLW  is  not  worthwhile  due  to  the  paucity  of  data,  and 
that  available  data  can  result  in  misleading  or  questionable  excursion 
rates.  Only  Wrightsville  and  Carolina  Beaches  have  high  temporal 
densities  of  data  points  for  each  MHW ,  MSL,  and  MLW  contour  and, 
consequently,  only  plots  from  these  beaches  were  redrawn  at  yet  a  larger 
scale  and  analyzed.  All  large-scale  plots  for  Wrightsville  Beach  and  a 
representative  set  from  Carolina  Beach  are  contained  in  Appendixes  A  and 
B,  respectively;  all  smaller  scaled  plots  for  Masonboro,  Kure,  and  Fort 
Fisher  Beaches  are  in  Appendixes  C,  D,  and  E,  respectively. 

Historic  events  which  may  have  affected  the  beach  erosion-accretion 
(excursion  distance)  are  indicated  on  each  excursion  distance  plot  for 
Wrightsville  and  Carolina  Beaches  (Figs.  16  and  17).  A  circle  is  placed 
on  a  data  point  measured  shortly  after  localized  storm  activity  (see 
Table  5),  and  an  arrow  is  placed  at  the  approximate  time  beach  fills 
were  completed.  The  same  profiles  in  Figures  11  (WB15)  and  12  (CB71) 
are  shown  in  Figures  16  and  17,  respectively,  drawn  at  the  larger  time 
scale  and  with  the  historic  events  indicated.  Excursion  rates  between 
the  beach  fills  (seasonally  averaged  response  shown  as  a  dashline)  can 
now  be  identified  and  quantified.  Localized  storms  account  for  many  of 
the  sudden  losses  in  beach  volume.  However,  some  erosion  (loss  of 
excursion  distance)  occurs  at  times  other  than  those  indicated  in 
Table  5,  possibly  due  to  localized  storms  of  lesser  magnitude,  but 
probably  due  to  erosion  from  swell  waves  generated  from  distant  storms. 

Sequential  beach  profiles  taken  between  January  1970  and  December 
1974  for  profile  WB15  are  presented  in  Appendix  F.  These  profiles  are 
presented  to  aid  the  reader  in  visualizing  the  postfill  response  of 
Wrightsville  Beach  and  thus  to  help  interpret  the  results  shown  in 
Figure  16. 

The  following  discussion  outlines  the  general  method  of  analysis 
used  on  all  excursion  distance  plots  for  Wrightsville  and  Carolina 
Beaches.  A  schematic  plot,  similar  to  the  MLW  excursion  distance  plot 
for  WB15  (Fig.  16),  is  used  as  an  example  and  is  shown  as  Figure  18. 
Section  IV. 4  contains  a  beach-by-beach  discussion  and  quantification 
detailing  the  effect  of  natural  and  manmade  influences  on  each. 

The  three  most  prominent  features  exhibited  by  Figures  16  to  18  are: 
(a)  the  long-term  erosional-accret ional  trend  is  approximately  constant 
(linear)  between  beach-fill  periods  with  minor  fluctuations  due  to 
seasonal  storm-induced  erosion  and  accretion  cycles;  (b)  the  placement 
of  a  fill  results  in  a  sudden  positive  spike  in  the  excursion  distances; 
and  (c)  immediately  following  a  significant  beach  fill,  loss  of  material 
occur 8  at  a  rapid  rate  which  gradually  decreases  to  equal  the  long-term 
recession  rate. 
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Not*:  Circles  Indicat*  profile*  measured  shorty  attar  a  local  storm 
Arrows  Indicats  tha  approximate  ttm*  at  which  beach  flits 
ware  placed. 


Figure  16.  Distance  from  the  base  line  to  stated  contours  at  WB  15. 


STANCE  IM1  DIS1RNCE  (M)  DISTANCE  (HI 


The  Long-term  change  for  most  beaches  in  the  study  is  negative, 
which  signifies  a  long-term  erosional  trend.  This  is  due  primarily  to 
the  inability  of  the  beach  to  return  to  its  original  position  after  a 
particularly  severe  winter  storm  period  or  after  a  very  severe  isolated 
storm  (e.g.,  a  hurricane  or  tropical  storm).  During  storm  activity, 
sediment  is  eroded  off  the  upper  section  of  the  beach  profile  and 
transported  either  alongshore  in  the  littoral  drift  or  offshore. 
Particularly  severe  storms  can  result  in  sediment  being  transported 
sufficiently  far  offshore  to  preclude  its  return  to  the  beach  face  under 
more  favorable  conditions,  thus  resulting  in  a  sediment  deficit  and, 
hence,  erosion.  Also  important  during  the  erosional  phase  of  beach 
behavior  is  the  continual  exposure  of  "fresh"  beach  sediment  which  may 
not  have  the  appropriate  sediment  distribution/characteristics  for  the 
dominant  wave  conditions.  This  means  that  under  erosional  conditions, 
sorting  losses  can  continually  occur  (resulting  in  long-term  losses), 
the  magnitude  of  which  is  dependent  upon  the  degree  of  mismatch  between 
the  distribution  of  the  exposed  sediment  to  that  which  is  more  suitable 
for  the  wave  conditions.  Another  cause  for  the  long-term  erosional 
problem  is  a  rise  in  sea  level  position.  Based  on  an  equilibrium  bottom 
profile,  Bruun  (1962)  quantified  the  volumetric  erosion  loss  per  unit 
Length  of  shoreline  (V)  as 


V  =  (e  +  d)  (X) 


(1) 


where  X  is  the  rate  of  shoreline  recession,  e  is  the  berm  crest  MSL,  and 
d  is  the  limiting  depth  between  nearshore  and  offshore  processes. 

Limiting  depth  (d)  is  approximately  -8.2  meters  (MSL)  based  on 
inspection  of  long  profiles  from  Wrightsville  and  Carolina  Beach  data. 
Horizontal  distances  to  this  depth  for  the  control  cells  are  presented 
in  Table  6.  The  rate  of  shoreline  recession  is  expressed  by 


ab 


(e  +  d) 


(2) 


where  a  is  the  rate  of  local  sea  level  rise,  and  b  is  the  distance  from 
the  initial  shoreline  to  the  limiting  depth. 


Table  6.  Volumetric  and  excursion  losses  due  to  rise  in  MSL. 


Littoral  Cell 

Distance  (b) 
to  1  im i t  i ng 
depth  of 
-8.2m 

Volumetric 
los  s/un  it 

Igth  of  beach 

(m-Vyr/m) 

Excursion  rate 
due  to  sea 
level  rise 
(m  /yr ) 

Wrightsville  Beach 

225 

0.83 

Masonboro  Beach 

210 

0.78 

Carolina  Beach 

190 

0.70 

Kure  Beach 

180 

0.67 

Fort  Fisher  Beach 

220 

0.81 

The  rise  in  MSL  during  the  study  period,  based  on  the  averaged 
trends  at  Portsmouth,  Virginia,  and  Charleston,  South  Carolina  (Hicks, 
1972),  was  approximately  0.37  centimeter  per  year.  The  computed  annual 
rate  of  volumetric  and  excursion  loss  due  to  the  rise  in  sea  level  for 
the  five  beaches  is  given  in  Table  6. 

The  rapid  loss  of  beach  material  immediately  after  the  placement  of 
a  beach  fill  can  be  split  into  two  components — a  long-term  component  due 
to  the  ongoing  long-term  processes,  and  an  initial  component  due  to 
enhanced  sorting  by  slope  readjustment.  The  continual  sorting  type 
losses  are  obviously  compounded  by  beach-fill  activity  when  sediment 
which  has  a  different  distribution  to  the  native  beach  sediment  is  used 
as  the  fill  material.  Not  only  is  the  magnitude  of  the  sorting  losses 
higher  because  of  the  generally  greater  mismatch  between  the  new 
distribution  and  the  desired  distribution,  but  also  the  rate  of  loss  is 
increased  due  to  the  increased  exposure  rate  to  wave  activity  as  a 
result  of  sediment  movement  due  to  slope  readjustment. 

The  long-term  component  can  be  represented  by  the  slope  of  the  line 
of  best  fit  through  all  data  points  after  time  t=t^  (Fig.  18),  such 
that  at  any  time,  t, 

lt  =  at  (3) 

where  lt  is  the  long-term  excursion  loss  (gain)  at  time  t,  and  a  is 
the  slope  of  the  linear  section  of  the  excursion  distance  plot. 

Data  from  this  study  indicated  that  after  1  to  2  years  following 
beach-fill  completion,  the  beach  face  generally  eroded  back  during  a 
winter  storm  period  to  its  approximate  prefill  position.  Both 
Figures  16  and  17  show  this  behavior  and  subsequent  accretion  of  the 
beach  face  during  the  ensuing  summer  period.  This  means  that  after 
approximately  2  years  most  of  the  beach-fill  material  has  been  exposed 
to  the  sorting  action  of  wave  activity  and  for  this  period  on  (i.e.,  the 
time  during  which  the  long-term  excursion  rates  were  calculated),  the 
enhanced  losses  due  to  the  sorting  of  beach-fill  material  should  have 
been  minimal. 

To  quantify  the  initial  loss  component,  the  long-term  component  was 
subtracted  from  the  excursion  distances  (shown  by  the  dashline  in 
Fig.  19).  The  time  scale  was  reset  to  zero  at  the  time  of  fill  (t=0), 
and  so  the  initial  loss  of  beach  fill  after  time  t  was  St .  Values 
of  S  (Fig.  19)  for  varying  time  increments  up  to  t=t^  were  plotted 
on  semilog  paper.  Figure  20  shows  the  results  of  these  plots  for  the 
MLW,  MSL,  and  MHW  excursion  curves  of  WB15.  The  results  from  this 
profile  are  typical  for  all  profiles  and  indicate  that  the  initial  loss 
component  due  to  sorting  and  beach-slope  adjustment  can  be  mathe¬ 
matically  represented  by  an  exponential  equation  of  the  form 

S  =  C  f.(l-10_kt)  for  0  <  t  <  t.  (4) 

t  l  t 
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where  k  is  the  slope  of  the  line  of  best  fit  of  the  semi  log  plot  of  S 
versus-  t,  f^  is  the  initial  fill  excursion,  £  is  the  fraction  of 
f t  lost  after  initial  losses  (i.e.,  at  t=t£),  and  St  is  the 
excursion  loss  at  time  t  due  to  sorting  and  slope  adjustments  of  a 
beach  fill. 

Note  that  the  exponential  form  of  equation  (4)  implies  that  the 
initial  losses,  although  very  small,  continue  indefinitely.  However, 
the  excursion  plots  indicate  that  after  1  to  2  years  the  excursion  loss 
sue  ic  slope  adjustment  and  initial  sorting  cannot  be  separated  from  the 
seasonal  and  long-term  losses.  Hence,  for  practical  reasons,  the 

. rial  loss  will  be  mathematically  considered  complete  when  95  percent 
of  s  fj,  is  lost  (i.e.,  at  t=ti). 

The  total  excursion  loss,  ,  at  time  t  after  fill  placement,  is 
the  sun  of  equations  (3)  and  (4). 

D  «  /  f.(l-10“kt)  +  at  (5) 

t  ^  l 

or,  the  total  beach  excursion  relative  to  the  prefill  position,  Et , 
at  ar.y  time  t  after  a  fill,  is 

Efc  =  f.  |"l-  C  +  <l(fktJ  -  at  (6) 

Equation  (6)  is  an  important  tool  which  can  be  used  to  evaluate 
historic  beach  fills  and  to  design  future  ones.  This  equation -can  be 
used  in  two  ways.  First,  if  a  given  design  lifetime  of  a  fill  is 
required,  substituting  Et=0  and  t  equal  to  the  desired  design  life, 
then  equation  (6)  is  solved  to  give  the  initial  fill  excursion  (and 
volume).  Second,  for  a  given  volume  of  fill,  or  alternatively,  for  a 
given  initial  excursion,  the  time  t=te  at  which  the  beach  returns  to 
its  prefill  position  (Et=0)  can  be  determined  (i.e.,  the  "useful 
life1'  of  the  fill  can  be  determined).  These  calculations  can  be  used  to 
quantify  the  effectiveness  and  value  of  a  given  beach  fill.  However, 
the  assumption  made  within  these  interpretations  of  equation  (6)  is  that 
the  beach  fill  has  lost  its  effectiveness  as  soon  as  the  beach  face 
between  the  MLW  to  MHW  contours  returns  to  its  initial,  prefill 
position.  It  must  be  noted  that  in  addition  to  providing  a  horizontal 
excursion  of  the  beach  face,  beach  fills  provide,  either  directly  or 
indirectly,  three  other  functions  which  retain  their  value  even  when  the 
initial  excursion  i3  lost.  The  direct  value  is  that  the  elevation  of 
the  berm(s)  and  sometimes  dunes  is  increased  during  beach-fill  opera¬ 
tions  so  that  a  larger  volume  of  material  seaward  of  the  backdune  is 
available  to  absorb  the  erosional  tendencies  of  storm  waves.  This  pro¬ 
vides  an  additional  degree  of  protection  to  the  backshore  which  was  not 
present  prior  to  the  fill  placement.  Indirectly,  beach  fills  result  in 
an  increase  in  sand  on  downdrift  beaches,  and  produce  slight  decreases 
in  the  nearshore  to  offshore  bathymetry  due  to  the  redistribution  of 
beach-fill  material  offshore  as  a  result  of  slope  readjustment  and 
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sorting.  These  decreased  depths  provide  an  added  measure  of  protect  on 
to  the  beach  by  forcing  waves  to  break  farther  offshore.  Individual 
designs  of,  and  the  nature  of  the  sediment  used  in  each  beach  fill, 
dictate  the  degree  to  which  these  factors  benefit  the  beach  area. 
Consequently,  they  will  not  be  further  addressed  in  this  analysis,  b-jt 
must  be  kept  in  mind  when  dealing  with  the  design  or  evaluation  of  a 
beach  fill. 

An  interesting  feature  of  Figure  20  is  the  relative  magnitude  of  :be 
k  values  (the  decay  rate)  of  the  MLW,  MSL,  and  MHW  curves.  The  greater 
the  k  value,  the  faster  the  rate  of  initial  loss  (erosion).  Conse¬ 
quently,  the  results  show  that  the  MHW  contour  eroded  at  a  faster  rate 
than  the  MSL  contour,  which  in  turn  eroded  at  a  faster  rate  than  the  MLW 
contour.  In  other  words,  the  slope  of  the  beach  face  readjusted  its’l: 
and  became  less  steep  during  the  initial  loss  period. 

4 .  Beach  Behavior  from  1965  to  1975. 

(a)  Wrightsville  Beach.  The  behavior  of  Wrightsville  Beach  in 
response  to  coastal  processes  during  the  1965  to  1975  decade  is  besr 
described  by  conveniently  dividing  Wrightsville  Beach  into  three 
sections — the  northern,  central,  and  southern  sections. 

The  northern  section  can  be  characterized  as  a  slowly  accreting 
beach  with  the  rate  of  accretion  falling  from  a  maximum  of  1.8  meter- 
per  year  at  Mason  Inlet  to  near  zero  about  1.75  kilometers  farther 
south.  Figure  21  shows  the  excursion  plots  for  WB3,  typical  of  the 
beach  behavior  in  this  northern  section.  Superimpose!  upon  the  aver 
accreting  excursion  is  a  seasonal  variation  of  appr oxin ately  20  met>  s 
The  minimum  excursion  distances  occur  during  the  first  three  (winter 
months  of  the  year  and  the  maximum  from  July  to  September.  Figure  2 
shows  that  the  beach  in  this  section  is  able  to  respond  to  storms, 
particularly  noted  are  those  in  February  and  March  of  1973,  and  to 
rebuild  itself  without  artificial  renourishment. 

Between  the  points  1.75  and  5  kilometers,  the  central  section  o 
Wrightsville  Beach  has  been  eroding  constantly  since  1965.  The  exc  r 
sion  plots  for  WB16  (Fig.  22)  are  typical  of  the  area  of  maximum  era;  on 
experienced  around  the  northern  area  of  the  town  of  Wrightsville  Beat  i. 
Beach  fills  in  1965,  1966,  and  1970  were  placed  to  protect  this  town 
however,  the  continued  high  erosion  rate  nullified  those  efforts.  T  o 
data  are  too  sparse  to  obtain  seasonal  variations  before  1970,  but  s  nee 
that  time  the  seasonal  excursion  within  the  central  section  was 
approximately  25  meters. 

The  behavior  of  the  southern  1 .5-kiloroeter  section  of  Wrightsvi  <• 
Beach  has  been  dominated  by  the  construction  of  the  northern  jetty  'o 
Masonboro  Inlet.  During  the  first  4  months  in  1966  (prior  to  the  19"6 
beach  fill),  the  nearshore  zone  of  the  beach  immediately  north  of  tii 
nearly  ^completed  jetty  accreted  by  up  to  40  meters,  especially  the 
and  MWL  contours  of  profiles  WB49  and  WB50.  This  accretion  fillet 
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CONTOUR  :  MHW 

Not*:  Circles  indicate  profiles  measured  shortly  after  a  local  storm. 

Arrows  indicate  the  approximate  time  at  which  beach  fills 
were  placed. 

Oashline  indicates  line  of  best  fit  of  average  excursion  distance. 


Figure  21 .  Distance  from  tha  bass  Ifna  to  stated  contours  at  WB  3 
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Notv:  Circles  indicate  profiles  measured  shortly  after  a  local  storm. 

Arrows  indicate  the  approximate  time  at  which  beach  fills 
were  placed. 

Oashline  indicates  line  of  best  fit  of  average  excursion  distance. 


Figure  22.  Distance  from  the  base  line  to  stated  contours  at  WB  16. 
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extended  northwards  with  time  into  the  area  of  beach  fill  and,  soon 
after  the  completion  of  the  jetty  in  spring  1966,  the  southern  end  of 
Wrightsville  Beach  had  accreted  by  approximately  30  to  40  meters.  From 
1 968  until  the  end  of  the  study  period,  the  accretion  fillet  underwent 
>nly  minor  changes  with  seasonal  fluctuations  of  15  to  25  meters. 

Figure  23  shows  typical  excursion  plots  of  WB47. 

The  long-term  excursion  rate  values  for  the  entire  beach  are  shown 
in  Table  7.  The  average  erosion  (excursion  loss)  per  year  along 
Wrightsville  Beach  due  to  the  rise  in  sea  level  is  0.10  meter  (see 
fable  6).  This  value  must  be  subtracted  from  the  measured  excursion 
rates  to  determine  the  average  annual  loss  of  beach  excursion  due 
primarily  to  longshore  processes.  These  values  are  shown  in  Table  7  and 
are  plotted  in  Figure  24. 

The  average  variation  in  seasonal  excursion  remained  fairly  constant 
ilong  the  entire  beach,  with  a  maximum  variation  occurring  at  MLW  and  a 
ainimum  at  MHW.  The  difference  in  the  seasonal  excursions  between 
1LW-MSL  and  MHW-MSL  gives  an  indication  of  the  average  change  in  beach- 
face  slope  from  winter  to  summer  beach  profiles.  Table  8  gives  the 
average  excursion  values  from  325  observations  along  Wrightsville  Beach; 
Figure  25  provides  a  visual  interpretation  of  the  relative  change  in 
seasonal  excursion  distances. 

There  were  insufficient  data  points  to  quantify  the  response  of 
Wrightsville  Beach  to  the  1965  and  1966  beach  fills.  However, 

Figure  26  shows  the  semilog  plots  of  the  initial  excursion  loss  after 
he  1970  beach  fill.  These  plots  show  the  combined  results  from  eight 
profiles  and  are  slightly  different  from  Figure  20.  The  values  of 
excursion  loss  at  time  t  after  beach-fill  placement  have  been  normalized 
t>y  dividing  them  by  the  total  initial  excursion  loss,  £f{,  and 
lence,  the  results  from  many  profiles  can  be  combined  to  compute  the 
average  exponential  decay  constant.  Table  9  gives  these  values  for  the 
1LW,  MSL,  and  MHW  contours,  together  with  values  of  C, ,  the  proportion  of 
he  MLW  to  MHW  fill  excursion  which  is  lost  due  to  sorting  and  slope 
idjustment,  the  initial  fill  excursion,  and  the  average  long-term  loss 
ate.  The  relative  differences  in  magnitude  of  the  k  values  for  the 
Lhree  contours  (shown  in  Table  9)  indicate  that  the  MSL  contour  eroded 
faster,  on  the  average,  than  either  the  MLW  or  MHW  contours,  thus 
producing,  as  expected,  a  concave  beach  profile.  The  average  long-term 
excursion  rate  of  -3.8  meters  (erosion)  per  year  for  all  three  contours 
indicates  that  once  long-term  slope  readjustments  occurred,  the  average 
beach  slope  did  not  change  from  year  to  year. 

(b)  Carolina  Beach.  Like  Wrightsville  Beach,  three  sections  of 
Jarolina  Beach  (northern  end,  north-central,  and  southern  half)  were 
affected  differently  by  the  action  of  the  coastal  processes  from  1965  to 
1975. 


The  northern  end  extends  from  Carolina  Beach  Inlet  southward  for 
1.5  kilometers  to  the  22-kilometer  point  (measured  from  the  northern 


TIME  IYERR3 


CONTOUR  :  MHW 

Noli:  Circles  indicate  profiles  measured  shortly  after  a  local  storm. 
Arrows  indicate  the  approximate  time  at  which  beach  fills 
were  placed. 

Oashline  indicates  line  of  best  fit  of  average  excursion  distance. 


Flour*  23.  Distance  from  the  baa*  tin*  to  stated  contours  at  WB  47, 
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Figure  24.  Comparison  of  measured  and  computed  volumetric 
change  along  Carolina  Beach. 


Table  8.  Seasonal  variation  in  MLW,  MSL,  and  MHW  position 
_ along  Wrightsville  Beach. 


Contour 

Avg  seasonal 

Excurs ion  minu s 

excur s ion 

MSL  excursion 

(m) 

(m) 

MLW 

28.9 

6.6 

MSL 

22.3 

0 

MHW 

20.6 

-1.7 

g 


23.9 


s  of  normalized  excursion  distance  versus  time  after  fill  placement 
LW,  and  MSL  contours  (1970  beach  fill). 


Figure  26.  Semi  log  plots  of  normalized  excursion  distance 
versus  time  after  fill  placement  for  MHW,  MLW, 
and  MSL  contours  (1970  beach  All).  -Continued 


Table  9.  Wrightsville  Beach,  1970  beach-fill  data. 


Contour 

Avg 

exponential 
decay  constant 
<k) 

Avg 

long-term 

excursion 

Initial 
beach-f il 1 
excursion  f- 
(m) 

t 

MLW 

0.58 

-3.8 

75.0 

0.80 

MSL 

0.63 

-3.8 

77.9 

0.79 

MHW 

0.59 

-3.8 

76.9 

0.82 
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limit  of  the  study  area).  Similar  to  the  northern  end  of  Wrightsville 
Beach,  this  section  of  Carolina  Beach  slowly  accreted  during  the  study 
period  with  a  maximum  rate  of  15  meters  per  year  at  the  tip  decreasing 
to  near  zero  at  22  kilometers.  As  shown  in  Figure  27,  this  area 
responded  naturally  to  storm-induced  erosion  and,  consequently,  no 
beach  fills  were  placed  during  the  study  period.  The  average  seasonal 
excursion  was  12.8  meters  for  the  northern  section. 

The  north-central  section  extends  from  the  22-  to  23 . 5-ki 1 ame t er 
points  and  encompasses  both  the  only  significant  change  in  beach 
orientation  along  Carolina  Beach  and  the  northern  end  of  the  town  of 
Carolina  Beach.  This  section  suffered  the  highest  measured  annual 
erosion  rate  of  the  entire  study  area,  and  estimates  of  that  rate  vary 
between  5  to  40  meters  per  year.  The  range  is  large,  and  errors  in  the 
estimation  of  the  excursion  rates  from  the  excursion  distance  plots 
probably  account  for  some  of  the  scatter  in  the  rate  values.  Because  of 
the  high  erosion  rates,  and  since  the  northern  end  of  the  town  of 
Carolina  Beach  is  exposed  to  this  erosion  (see  Fig.  4),  six  beach  fills 
were  placed  in  this  section  between  1965  and  1971,  three  of  which  were 
connected  with  the  experimental  deposition  basin  in  the  throat  of 
Carolina  Beach  Inlet.  The  excursion  distance  plots  for  CB64  (Fig.  28) 
reveal  rapid  erosion  after  each  beach  fill  and  the  continued  loss  of 
beach  material  despite  the  beach-fill  activities.  The  seasonal 
excursion  distance  within  this  area  is  about  19.5  meters. 

The  southern  half  of  Carolina  Beach  experienced  mild  erosion  rates 
of  approximately  5  meters  per  year.  Beach  fills  in  1965  and  1971 
provided  protection  to  the  southern  end  of  the  Carolina  Beach  township 
because  the  net  excursion  in  1974  was  still  positive;  i.e.,  more  sand 
was  placed  on  the  beach  by  the  beach-fill  projects  than  was  eroded  away 
during  the  1965-74  period.  Figure  29  shows  an  example  of  the  excess  in 
excursion  distance  for  CB119  and  also  shows  that  the  average  seasonal 
variation  along  this  section  is  relatively  small  with  a  mean  value  of 
approximately  7.6  meters. 

The  long-term  excursion  rates  for  the  entire  beach  are  shown  in 
Table  10.  The  representative  value  of  average  annual  excursion  loss 
along  Carolina  Beach  due  to  the  rise  in  sea  level  is  0.09  meter  (see 
Table  6).  This  value  must  be  subtracted  from  the  measured  excursion 
rates  to  determine  the  annual  excursion  loss  due  to  longshore  processes. 
Representative  values  are  given  in  Table  10,  and  a  complete  set  along 
Carolina  Beach  is  plotted  in  Figure  30. 

Table  11  shows  the  average  MLW,  MSL,  and  MHW  seasonal  excursion 
values  for  the  entire  beach  and  the  relative  differences  in  seasonal 
variation  between  these  contours.  The  average  change  in  beach  slope  at 
MSL  from  a  summer  profile  to  a  winter  profile  was  0.2°,  i.e.,  1  on  286. 
Figures  31  and  32  show  the  semilog  plots  of  the  normalized  initial 
excursion  loss  values  versus  time  after  fill  placement  for  the  1965  and 
1971  beach  fills,  respectively.  Since  there  is  a  lack  of  data  for  the 
1971  fill,  all  MLW,  MSL,  and  MHW  values  from  profile  CB93  were  combined 
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TIME  (T  ERR) 

CONTOUR  :  MHW 

Noli:  Circles  indicate  profiles  measured  shortly  after  a  local  storm. 

Arrows  indicate  the  approximate  time  at  which  beach  fills 
were  placed. 

Dashline  indicates  line  of  best  fit  of  average  excursion  distance. 


Figure  27.  Distance  from  the  base  line  to  stated  contours  at  CB  2 
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CONTOUR  :  MHW 

Note:  Circles  indicate  profiles  measured  shortly  after  a  local  storm. 
Arrows  indicate  the  approximate  time  at  which  beach  fills 
were  placed. 

Dashline  indicates  line  of  best  fit  of  average  excursion  distance. 

Flyura  29.  Distance  from  tho  basa  lino  to  stated  contours  at  CB  119. 
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ANNUAL  VOLUME  CHANGE  PER  UNIT  LENGTH  OF  BEACH 


Table  10.  Average  long-term  excursion  rates 


along  Carolina  Beach. 


Profile 

station 

Distance  from 
north  study 
boundary 
(km) 

Avg 

excursion 
rate 
(m/yr ) 

Avg  excursion 
rate  due  to  long¬ 
shore  processes! 
(m/yr  ) 

CBl2 

20.64 

9.0 

9.1 

CB2 

20.94 

15.0 

15.1 

CB10 

21.70 

6.0 

6.  1 

CE15 

22.15 

-12.9 

-12.8 

CB  16 

22.23 

-17.9 

-17.8 

CB21 

22.35 

-20.9 

-20.8 

CB32 

22.48 

-22.5 

-22.4 

CB40 

22.53 

-10.4 

10.3 

CB44 

22.54 

-5.3 

-5.2 

CB53 

22.59 

-18.5 

-18.4 

CB6I 

22.72 

-5.8 

-5.7 

CB71 

23.03 

-27.4 

-27.3 

CB84 

23.36 

-19.2 

-19.1 

CB93 

23.64 

-9.2 

-9.1 

CB96 

24.24 

-5.8 

-5.7 

CB99 

25.20 

-4.9 

-4.8 

CB106 

26.07 

-3.4 

-3.3 

CB  1 1  7 

26.56 

-1 .0 

-0.9 

CB  1 19 

26.68 

-0.8 

-0.7 

^Based  on  a  0 .09-meter  loss  in  excursion  due  to  a  rise  in  sea  level. 
^Profiles  within  inlet  shoals. 
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Figure  30.  Comparison  of  measured  and  computsd  volumatrlc 
changa  along  Wrlghtsvllla  Baach. 


NORMALIZED  EXCURSION  DISTANCE 


Table  11.  Seasonal  variation  in  MLW,  MSL,  and  MHW 
positions,  Carolina  Beach. 


Contour 

Avg 

seasonal 

excursion 

(tn) 

Excursion  minus 

MSL  excursion 

(m) 

MLW 

18.2 

0.4 

MSL 

17.8 

0 

MHW 

16. A 

-1.4 

Figure  31.  Semllog  plots  of  normalized  excursion  distance  versus  time  after 
fill  placement  for  1965  beech  fill. 
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Figure  32.  Semilog  plots  of  normalizod  excursion  dlstsnce  versus 
time  sfter  fill  plscement  for  1971  beech  fill. 
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to  calculate  the  exponential  decay  constant  for  the  sorting  and  slope 
adjustment  losses.  Table  12  contains  all  relevant  data  for  the  1965  and 
1971  beach  fills  that  could  be  confidently  extracted  from  the  excursion 
distance  plots. 


Table  12.  1965  to  1971  beach-fill  data. 


Beach  fill 


Carolina  Beach. 


Avg  exponential 

decay  count  (k) 

MLW 

MSL 

MHW 

Avg 

0.83 

0.98 

0.70 

0.84 

n  o  *1 

Avg 

initial 

fill 

excursion 

(m) 


Avg 

long-term 
excurs ion 

(tn/yr ) 


(c)  Masonboro,  Kure,  and  Fort  Fisher  Beaches.  Because  of  insuffi¬ 
cient  and  nonconsistent  temporal  distribution  of  excursion  distance 
data,  beach  response  in  terms  of  long-term  eros ional-accretional  rates, 
beach  fills,  and  storm  events  cannot  be  described  for  Masonboro,  Kure, 
or  Fort  Fisher  Beaches.  Therefore,  only  a  brief  statement  concerning 
the  relative  difference  in  excursion  distance  between  the  first  and 
final  data  points  can  be  made;  however,  because  of  seasonal  variation 
=»od  possible  poststorm  excursions,  even  this  may  be  misleading. 


From  1966  to  1973,  the  erosional  loss  at  Masonboro  Beach  was 
generally  10  to  30  meters.  However,  two  profiles  (MB2  and  MB5) ,  which 
are  located  in  the  vicinity  of  the  only  significant  change  in  beach 
angle  along  Masonboro  Beach,  show  losses  of  80  to  100  meters.  The 
excursion  differences  for  most  profiles  fall  within  the  possible  range 
of  seasonal  or  poststorm  excursion  ranges  and,  consequently,  the  actual 
long-term  loss  on  Masonboro  Beach  may  not  be  reflected  by  the  above 
values . 


The  availability  of  excursion  distance  data  for  Kure  Beach  and  Fort 
Fisher  Beach  is  even  less  than  that  for  Masonboro  Beach,  with  data 
collected  only  frcxn  late  1969  to  early  1973.  Differences  in  excursion 
positions  between  those  dates  for  both  beaches  vary  from  +5  to 
-20  meters,  but  again,  estimated  seasonal  variation  from  two  profiles  of 
10  to  15  meters  makes  any  conclusion  on  the  long-term  response  of  these 
beaches  impossible. 


V.  LONGSHORE  SEDIMENT  TRANSPORT  ANALYSIS 


1 .  Introduction. 

The  procedure  to  mathematically  predict  the  volume  of  sediment  in 
the  littoral  drift  requires  knowledge  of  the  magnitude  and  direction  of 
the  energy  flux  due  to  waves  breaking  along  the  study  area  beaches.  To 
determine  this  quantity,  a  wave  climate  representative  of  the  annual 
wave  conditions  measured  or  experienced  in  offshore  waters  must  be 
established.  The  wave  climate,  in  this  case  in  the  form  of  a  set  of 
wave  heights  with  different  periods  and  directions,  must  be  "routed" 
towards  shore  by  a  wave  refraction  model  until  the  waves  break  on  or 
near  the  beach.  Information  on  their  breaking  angles  (relative  to  the 
beach  orientation),  breaking  wave  heights,  and  wave  speed  at  breaking 
are  determined  and  used  to  establish  the  longshore  components  of  the 
energy  flux  for  both  the  northerly  and  southerly  directions. 

The  quantity  of  sediment  carried  by  the  littoral  drift  in  each 
direction  is  found  by  multiplying  the  magnitude  of  the  energy  flux  by  a 
conversion  factor  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  1977).  However,  uncertainty  exists  in  the  exact  value 
of  that  factor  (Vitale,  1980),  and  therefore,  it  will  be  recalculated 
for  this  study  area  by  comparing  the  known  time  rate  of  volumetric 
change  at  Wrightsville  Beach  and  Carolina  Beach  to  the  predicted  values 
of  the  energy  flux  at  those  beaches.  The  recomputed  conversion  factors 
will  be  used  to  estimate  the  annual  northerly  and  southerly  longshore 
transport  quantities  and  the  volume  of  material  lost  into  the  adjacent 
inlets. 

2.  Wave  Refraction  Analysis. 

(a)  Wave  Climate.  Wave  climate  was  determined  from  a  joint 
probability  evaluation  of  wave  gage  data  at  Johnnie  Mercer's  Pier  and 
wave  observation  data  from  Wrightsville  Beach.  The  directional 
distribution  of  wave  height  and  wave  period,  calculated  from  the  wave 
observation  data,  was  assumed  to  hold  for  the  Johnnie  Mercer's  Pier 
data.  Consequently,  wave  angles  at  the  gage  were  statistically 
correlated  to  the  wave  observation  data  observations.  The  SSMO  and 
Frying  Pan  Shoals  wave  data  were  not  used  due  to  a  lack  of  confidence  in 
data  recording  (Harris,  1972). 

Under  random  sea  conditions,  the  distribution  of  the  values  for  wave 
height,  period,  and  direction  is  continuous.  However,  to  perform  the 
wave  refraction  analysis,  a  representative  set  of  wave  height,  period, 
and  direction  conditions  was  needed.  Consequently ,  the  distribution  of 
wave  height  was  divided  into  three  ranges  and  the  period  into  six  groups 
with  midrange  values  of  3,  6.3,  8.3,  10.3,  12.3,  and  16  seconds.  The 
angles  of  wave  approach  were  also  divided  into  four  sectors  (northeast, 
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east,  southeast,  and  south),  with  the  wave  statistics  from  the  inter¬ 
mediate  directions  (north-northeast,  east-northeast,  etc.)  being  incor¬ 
porated  proportionately  into  the  four  primary  directions.  Figure  33 
shows  these  approach  angles  relative  to  the  shoreline  orientation. 

The  distribution  of  wave  height  was  converted  to  an  equivalent 
distribution  of  wave  energy  (wave  height  squared)  and  divided  into  three 
ranges.  The  wave  height  corresponding  to  each  of  the  midrange  values  of 
wave  energy  was  then  determined.  The  offshore  wave  height  and  approach 
angle  corresponding  to  each  of  the  three  nearshore  wave  heights  were 
calculated  for  each  period  and  nearshore  angle  condition.  Both  the 
offshore  wave  direction  and  refraction  coefficients  were  determined  by 
using  Snell's  Law,  and  the  shoaling  coefficients  were  calculated  by  the 
ratio  of  nearshore  and  offshore  depths.  The  offshore  wave  heights  cor¬ 
responding  to  each  of  the  three  nearshore  wave  heights  were  calculated 
by  dividing  the  nearshore  height  by  the  product  of  the  refraction, 
shoaling,  and  friction  coefficients.  Explanation  of  the  development  of 
the  friction  coefficient  is  detailed  later  in  Section  (c).  Ibe  three 
offshore  wave  heights  used  in  the  analysis  were  0.32,  1.40,  and 
2.47  meters. 

The  probability  of  occurrence  (expressed  as  a  percentage)  of  a  wave 
approaching  the  study  area  from  each  of  the  four  directions,  with  a  wave 
height  and  period  falling  within  one  of  the  three  height  ranges  and  six 
period  ranges  (i.e.,  72  different  cases),  was  calculated  from  the  data 
sets  for  each  season;  i.e.,  winter  (December,  January,  and  February), 
spring  (March,  April,  and  May),  summer  (June,  July,  and  August),  and 
fall  (September,  October,  and  November).  This  information  is  presented 
in  Table  13. 

The  percentage  of  occurrence  of  many  of  the  wave  heigh t-per iod- 
direction  combinations  is  less  than  one.  To  reduce  excessive  and 
unnecessary  analysis  costs,  it  was  decided  that  satisfactory  results 
could  be  achieved  by  using  only  enough  wave  combinations  so  that,  for 
each  season,  95  percent  of  occurrence  by  wave  energy  of  all  possible 
combinations  of  height,  period,  and  direction  was  modeled.  Selection  of 
seasonal  wave  types  was  based  on  the  summation  of  percentage  of 
occurrence  by  wave  energy  of  those  wave  conditions  with  the  highest 
percentage  until  the  95-percent  criterion  was  satisfied.  Summation  to 
95  percent  by  wave  energy  resulted  in  a  representation  of  the  wave 
climate  by  approximately  98  percent  of  the  observed  wave  types.  Table  14 
shows  the  offshore  wave  climate  chosen  to  represent  the  average  seasonal 
conditions  measured  along  the  study  area.  The  average  annual  climate  is 
represented  by  the  arithmetic  average  of  the  seasonal  values  for  each 
combination  of  wave  height,  period,  and  direction. 

The  final  step  in  the  selection  of  the  wave  climate  data  was  a 
calibration  check  using  the  wave  refraction  model.  The  annual  wave 
climate  sets  were  refracted  toward  shore  and  combined  according  to  their 
percentage  of  occurrence  (see  Section  V,  3).  The  directional 
distribution  of  the  wave  energy  at  Wrightsville  Beach  was  compared  to 


il  percent  of  occurrence  for  that  specific  wave  condition 


the  measured  distribution  calculated  from  the  wave  observation  data. 
Considering  the  errors  inherent  in  tin  visual  data  collection  method,  in 
the  data  analyses  techniques,  and  errors  resulting  from  presenting  the 
continuous  distribution  of  wave  approach  angles  as  approach  sectors. 
Table  15  shows  a  favorable  comparison. 


Table  15.  Predicted  and  measured  distribution  of  wave  energy  at 
Wrightsville  Beach. 


(b)  Bathymetric  Data.  The  wave  refraction  model  requires  knowledge 
of  the  general  bathymetry  offshore  from  the  study  area  to  accurately 
refract  the  approaching  wave  sets.  The  bathymetric  data  was  provided  on 
a  150-meter  (500-foot)  square-grid  spacing  which  extended  from  the  MLW 
position  of  the  shoreline  to  a  depth  of  appr oximat ley  20  meters 

(65  feet),  15  kilometers  (9.4  miles)  offshore.  The  nearshore  depths 
were  interpolated  from  the  long  beach  profiles  and  the  greater  offshore 
depths  were  measured  from  1978  National  Ocean  Survey  (NOS)  nautical 
charts . 

The  offshore  bathymetry  of  the  study  area  is  quite  irregular  and  a 
qualitative  graphical  representation  of  it  is  shown  in  Figure  34.  This 
figure  is  a  three-dimensional  line  drawing  display  of  the  data  generated 
by  a  computer  graphics  program,  and  consequently  the  offshore 
representation  is  quite  accurate.  However,  the  interpolation  scheme 
used  by  this  program  distorted  the  shoreline  position,  and  a  dot  screen 
pattern  has  been  included  to  alleviate  this  visual  distraction. 

(c)  Wave  Refraction  Model.  The  numerical  model  used  for  the  wave 
refraction  analysis  is  a  modified  version  of  the  wave  refraction  model 
developed  by  Dobson  (1967).  Dobson's  model  requires  the  wave  ray  to 
originate  in  deep  water,  a  condition  which  is  not  always  practical  (or 
economical  relative  to  computer  costs)  for  long-period  waves.  There¬ 
fore,  a  subroutine  was  added  to  account  for  the  refraction  and  shoaling 
of  the  wave  ray  which  occurs  in  the  deeper  offshore  regions.  This 
routine  assumes  that  bathymetry  in  the  offshore  region  has  straight 
and  parallel  contours.  Snell's  law  is  used  to  compute  the  refraction 
coefficient  and  the  change  in  the  wave  angle  at  an  economically  more 
reasonable  "offshore”  boundary  for  the  model.  The  partially  refracted 
wave  ray  is  then  used  as  the  starting  condition  for  Dobson’s  numerical 
model  which  integrates  the  wave  ray  through  shallower  regions  toward  the 


MSL  shoreline.  For  this  study,  the  numerical  model  offshore  boundary 
extended  to  about  the  20-meter  (65-foot)  depth  contour  (MSL),  about 
15  kilometers  (9.4  miles)  offshore. 

A  second  modification  to  the  original  program  was  the  addition  of  a 
subroutine  to  account  for  energy  losses  due  to  friction.  The  wave 
height,  H,  at  any  point  along  the  wave  ray  can  be  represented  by 

H  =  H0-Kr-Ks-Kf  (7) 

where  HQ  is  the  deepwater  wave  height,  Kr  is  the  refraction 
coefficient,  Kg  is  the  shoaling  coefficient,  and  K£  is  the  friction 
coefficient . 

Dobson's  (1967)  original  model  calculated  both  the  refraction  and 
shoaling  coefficients.  The  additional  subroutine  calculates  the  fric¬ 
tion  coefficient  by  integrating  an  expression  developed  by  Skovgaard, 
Jonsson,  and  Bertelson  (1975)  along  the  wave  ray  from  deep  water  to  the 
point  of  interest  (optionally  the  point  of  wave  breaking).  The  integra¬ 
tion  is  carried  out  using  a  trapezoidal  integration  scheme.  The  local 
bottom  friction  factor  is  calculated  from  the  local  wave  conditions  by  a 
numerical  algorithm  developed  by  Fritsch,  Shafer,  and  Crowley  (1973). 

The  expression  for  the  wave  friction  coefficient,  as  given  by  Skovgaard, 
Jonsson,  and  Bertelsen,  further  requires  a  value  for  the  equivalent 
(Nikuradse)  bottom  roughness.  A  field  observation  on  a  sandy  coast  by 
Iwagaki  and  Kakinuma  (1963)  found  that  the  bottom  roughness  ranged  from 
l  to  2  centimeters.  For  this  study,  the  value  of  equivalent  bottom 
roughness  was  determined  from  the  calibration  of  offshore  SSMO  wave 
height  (wave  energy)  data  which  had  been  routed  inshore  to  wave  height 
(wave  energy)  data  measured  at  Johnnie  Mercer's  Pier  gage.  Although 
some  uncertainty  exists  with  the  SSMO  data,  as  noted  in  Section  2(a),  it 
was  used  here  in  a  simple  test  to  determine  whether  or  not  the 
literature  values  for  bottom  roughness  were  applicable  on  this  part  of 
the  coast.  A  value  of  1.5  centimeters  gave  the  best  results  for  the 
comparison  of  computed  and  measured  wave  energy  at  the  beach,  and  this 
value  falls  within  Iwagaki  and  Kakinuma's  range  of  values. 

The  effect  of  including  bottom  friction  in  the  wave  refraction  model 
is  a  reduction  in  the  wave  height  and,  therefore,  wave  energy  as  the 
wave  ray  progresses  into  shallow  water.  It  has  no  effect,  within  the 
limits  of  the  linear  theory  used  by  Dobson  (1967),  on  the  direction  of 
wave  propagation;  however,  reduction  of  the  wave  height  does  affect 
breaking  conditions,  as  a  wave  with  a  reduced  height  can  propagate 
closer  to  shore  before  breaking.  For  waves  in  shallow  water,  solitary 
wave  theory  defines  the  breaking  condition 


5-  =  0.78  (8) 

a 

where  H  is  the  local  wave  height,  and  d  is  the  local  water  depth. 

The  third  modification  to  Dobson's  model  was  a  routine  to  stop 
integration  of  the  wave  ray  when  the  ratio  of  wave  height  to  local  water 
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depth  exceeds  0.78.  To  determine  the  depth  at  any  point  along  the  wave 
ray,  the  model  uses  an  algorithm  which  fits  a  polynomial  to  the  depth  of 
the  surrounding  square  of  eight  grid  points  (relative  to  that  wave  ray). 
Under  the  rapidly  varying  bathymetric  conditions  which  exist  within  the 
study  area,  the  algorithm  often  computed  nonrepresentative  depth  values 
which  in  turn  resulted  in  offshore  wave  breaking  and  caustic  (wave 
crossing)  conditions.  To  help  alleviate  this  problem,  the  depth  grid 
spacing  was  increased  from  150  meters  (500  feet)  to  300  meters 
(1,000  feet),  and  this  modification  resulted  in  a  significant  reduction 
in  the  number  of  offshore  caustics  and  wave  breaking.  In  addition  to 
this  problem,  diffraction  (i.e.,  the  lateral  spreading  of  energy  along 
the  crest  of  a  wave),  an  important  process  in  "smoothing-out"  peaks  in 
wave  energy  (and  height),  is  ignored  by  Dobson's  model. 

Figures  35  and  36  are  two  computer-generated  wave  refraction 
diagrams  for  a  wave  approaching  from  the  east  with  an  offshore  wave 
height  of  1.4  meters  and  a  period  of  10.5  seconds.  Figure  35  shows  that 
many  of  the  wave  rays  cross  before  reaching  the  beach  or  break  offshore. 
Since  each  wave  ray  is  propagated  independently  toward  the  shoreline, 
the  model  is  "unaware"  of  the  possibility  that  any  two  or  more  wave  rays 
may  cross.  Linear  wave  theory  is  not  valid  under  these  conditions; 
therefore,  all  wave  rays  which  crossed  before  reaching  breaking  condi¬ 
tion  must  be  eliminated  from  the  analysis.  Figure  36  shows  the  same 
wave  propagation  as  in  Figure  35;  however,  all  crossed  wave  rays  have 
been  eliminated.  The  energy,  and  therefore,  wave  properties  like 
height,  celerity,  and  angle  along  a  wave  crest  between  two  adjacent 
noncaustic  rays,  was  assumed  to  be  proportional  to  the  energy  values  of 
these  noncaustic  rays.  Hence,  breaking  wave  conditions  at  all  locations 
along  the  beach  were  found  by  linearly  interpolating  the  values  between 
adjacent  noncaustic  wave-ray  locations. 

Another  shortcoming  of  Dobson's  (1967)  model  is  that  the  influence 
of  tidal  jets  and  currents  near  inlets  on  wave  refraction  is  not 
considered.  Together  with  the  fact  that  bathymetric  changes  are  rapid 
in  the  vicinity  of  inlets,  the  resulting  values  of  wave  height,  angle, 
and  celerity  at  those  locations  must  be  considered  with  some  skepticism. 

Computer  plots  showing  the  results  of  the  refraction  analysis  for 
1.4-meter  waves  for  each  wave  period  and  for  all  four  wave  approach 
angles  are  contained  in  Appendix  G.  The  difference  between  the  results 
of  waves  having  the  same  period  and  approach  direction,  but  differing  in 
height,  is  simply  a  slight  difference  in  the  breaking  position  of  the 
wave  along  the  same  wave-ray  path. 

3 .  Energy  Flux  Computation. 

The  longshore  component  of  wave  energy  flux,  ,  is  defined  as 
(U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center, 
1977;  Vitale,  1980) 


?! 


sin  2  n 


(9) 
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crossed  wave  waves  eliminated. 


where  H  is  the  wave  height,  Cg  is  the  wave  group  velocity,  and  a  is  the 
angle  the  wave  crest  makes  with  the  shoreline.  Usually  the  breaking 
wave  characteristics  (Hj>,  Cg^,  and  a are  used  to  represent  the  wave 
energy  flux  entering  the  surf  zone. 

Each  wave  type  was  refracted  toward  shore  by  the  refraction  model. 
The  breaking  wave  values  of  H^,  Cg,  and  approach  angle,  were  deter¬ 
mined  at  each  breaking  wave-ray  location,  and  then  interpolated  at  beach 
stations  every  250  meters  along  the  study  area.  The  shoreline  (plan) 
angle  at  each  of  these  250-meter  locations  was  measured  from  aerial 
photos  and  the  value  of  a  then  determined.  The  longshore  component  of 
wave  energy  flux  at  breaking  was  calculated  using  equation  (9)  at  each 
250-meter  beach  station,  and  was  then  multiplied  by  that  wave  type's 
percent  occurrence.  A  positive  value  of  Pj^  represented  a  component 
of  wave  energy  flux  in  a  southerly  direction  and  a  negative  value 
represented  a  component  in  the  northerly  direction. 

As  each  wave  type  was  refracted  toward  shore,  and  the  longshore 
component  of  wave  energy  flux  was  calculated,  the  percent  contribution 
to  either  the  northerly  or  southerly  components  of  the  annual  longshore 
flux  was  summed,  by  direction,  with  the  contribution  from  the  other  wa *e 
types.  The  resulting  totals  at  each  250-meter  beach  station  represent 
the  northerly  and  southerly  longshore  components  of  the  annual  wave 
energy  flux. 

The  spatial  variation  of  these  totals  was  significant,  and  the 
sudden  changes  in  magnitude  were  not  representative  of  the  actual  energy 
flux  conditions.  Several  factors  which  contributed  to  this  problem 
were: 

(a)  The  refraction  model  used  a  static  representation  of  shoreline 
conditions  and  bathymetry.  As  soon  as  a  concentration  of  wave 
energy  in  shallow  water  occurs  in  the  prototype,  erosion 
results  and  bathymetry  changes  to  reduce  the  energy  concen¬ 
tration;  i.e.,  nature  tends  to  smooth  out  sudden  changes  in 
concentrations  of  wave  energy,  but  the  model  cannot. 

(b)  The  resolution  of  the  computational  grid  cells  close  to  the 
beach  were  not  fine  enough  to  allow  for  the  rapid  changes  in 
bathymetry  and  beach  planform. 

(c)  The  energy  flux  values  are  proportional  to  the  product  of  the 
sine  and  cosine  values  of  the  wave  approach  angle  relative  to 
the  beach  shoreline.  Consequently,  subtle  errors  in  offshore 
angles  can  result  in  significant  errors  in  the  energy  flux 
computation  at  the  beach  face. 

(d)  Diffraction  effects  and  the  influence  of  tidal  currents  were 
not  included. 
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To  overcome  these  problems,  i.e.,  to  remove  the  rapid  fluctuations 
without  significantly  altering  the  longer  terra  treads,  a  nine-point 
running  filter  was  applied  to  the  results  of  the  energy  flux  computa¬ 
tions.  The  running  filter  averages  the  values  from  nine  points  (in  this 
case,  nine  250-meter  points  are  equivalent  to  averaging  over  a 
2-kilometer  stretch  of  beach)  and  assigns  that  average  to  the  middle 
point.  The  filter  is  then  moved  to  the  next  (middle)  point  and  averages 
its  value  with  the  four  values  on  either  side,  etc. 

Figure  37  shows  the  filtered  results  of  the  northerly  and  southerly 
components  of  the  annual  longshore  energy  flux;  Figure  38  combines  both 
components  and  shows  the  net  annual  longshore  energy  flux  acting  along 
the  study  area. 

4.  Longshore  Sediment  Transport  Model. 

The  accepted  practice  for  computing  the  longshore  sediment  transport 
rate  has  been  to  use  an  empirical  relationship  between  the  longshore 
component  of  the  energy  flux  entering  the  surf  zone  and  the  volume  of 
sand  moved.  This  dimensional  relationship  is  given  in  the  Shore  Pro¬ 
tection  Manual  (SPM)  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  1977)  and  can  be  expressed  as 

Q  —  =  1,288  PX 

yr  N-yr  xs  s-M 

or 

Q  ^  •=  7,500  ^  Pl 

yr  lb-yr  As  ft-s 

where  P^  is  the  energy  flux  factor  and  0  is  the  longshore  sediment 
transport?  rate.  This  equation  was  developed  from  field  observations  in 
which  wave  height  characteristics  were  represented  by  only  one  value — 
the  significant  wave  height. 

In  this  study,  actual  longshore  energy  flux  components  were  cal¬ 
culated  for  a  set  of  wave  types  which  were  subsequently  summed  together 
according  to  their  percent  occurrence.  Consequently,  this  calculation 
of  the  longshore  energy  flux  is  not  compatible  with  equation  (10)  above; 
hence,  the  dimensional  constants  given  in  the  SPM  cannot  be  directly 
applied  or  compared.  Jarrett  (1977)  performed  a  refraction  analysis 
similar  to  that  performed  in  this  study  and  found  a  value  for  the 
constant  by  correlating  measured  volumetric  changes  along  Wrightsville 
Beach  to  computed  energy  flux  values  at  each  end  of  the  beach. 

Jarrett's  successful  results  showed  that  the  sane  type  of  relationship 
which  is  given  in  the  SPM  exists  between  the  computed  values  of  the 
longshore  energy  flux  and  the  sediment  transport  rates.  Therefore,  that 
relationship  is  used  in  this  study  and  is  expressed  by 


(10a) 


(10b) 
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Figure  38.  Net  annual  longshore  energy  flux  along  study 


where  a  is  the  number  of  wave  types  used  to  represent  the  seasonal  or 
annual  wave  climate,  P^.  is  the  longshore  component  of  wave 
energy  flux  (at  breaking),  p^  is  the  percent  occurrence  of  that  wave 
type,  Q  is  the  long-term  longshore  sediment  transport  rate,  and  )3  is 
the  dimensional  constant  (found  from  correlation)  relating  Q  to  Pk. 

The  dimensions  of  each  term  are  shown  in  brackets. 

A  sediment  budget  approach  can  be  used  for  the  correlation  of  Q  and 
Pi*  For  a  beach  cell,  as  shown  below,  Qin  represents  all  long-term 
sources  of  sediment  supplied  into  the  cell-per-unit  time  and  Qout 
all  long-term  losses  from  the  cell-per-unit  time.  The  difference, 
Qout~Qin«  represents  the  long-term  change  in  beach  volume  for 
'that  cell. 

The  longshore  components  of  wave  energy  flux,  as  calculated  in 
Section  V,3,  are  Pj,  and  P^_  and  their  respective  beach  coordinates  are 
Xi  and  X2. 


X1  , _ x2  v 

Qin - ^  1  beach  cell] - >  Qout 

pll  Pl2 

From  equation  (11),  QoUt-Qin  -  /9  [pli~pl2]-  Let  q/L  be  the  long 
-term  erosion  or  accretion  rate  per  unit  length  of  beach,  then 
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At  any  point  along  the  beach,  /9  can  be  determined  from  the  ratio  of  the 
long-term  erosion-accretion  rate  to  the  spatial  gradient  of  the 
longshore  component  of  wave  energy  flux. 


Values  of  measured  q^  were  taken  from  all  profiles  along  the 
beaches  away  from  the  immediate  area  of  inlet  influence.  Unfortunately 
due  to  the  insufficient  temporal  and  spatial  distribution  of  profile 
data,  volumetric  change  data  for  Masonboro,  Kure,  and  Fort  Fisher 
Beaches  were  not  calculated.  Only  values  for  Wrightsville  and  Carolina 
Beaches,  in  Tables  7  and  10,  were  compared  to  predicted  values.  A  plot 
of  q^  and  /3  (dPji/dX),  versus  beach  distance  X,  was  drawn  by 
choosing  a  value  of  /3  which  produced  the  best  correlation  between  the 
two  lines.  To  eliminate  sudden  computational  fluctuations  before 
comparison  with  measured  q^  values,  the  (3  dP^/dX  values  were 
filtered  to  produce  smoothly  varying  distribution. 

Figures  39  and  40  show  the  results  of  these  comparisons  for 
Wrightsville  Beach  and  Carolina  Beach,  respectively.  Although  consid¬ 
erable  scatter  in  the  values  of  q^  is  obvious,  especially  along  the 
northern  Carolina  Beach  region,  the  general  trends  of  both  the  computed 
and  measured  volumetric  change  values  are  similar  along  each  beach. 

Within  the  limitations  of  the  analysis,  it  appears  that  a  value  of 
/3  =300  m^-s/N-yr  provides  the  best  fit  for  Wrightsville  Beach  with  a 
data  scatter  of  +33  percent.  For  Carolina  Beach,  the  best-fit  value  is 
/3  =900  m^-s/N-yr  with  a  data  scatter  of  +66  percent.  These  results 
as  summarized  in  Table  16,  show  a  large  possible  range  in  values  of  /3 . 
Assuming  that  equation  (11)  is  a  valid  representation  of  the  relation¬ 
ship  between  the  longshore  sediment  transport  rate  and  the  longshore 
component  of  wave  energy  flux,  then  two  possible  conclusions  can  be 
made.  First,  the  value  of  $  is  highly  localized  and  strongly  dependent 
on  the  local  physical  characteristics  of  the  beach  and  sediment 
properties.  Table  3  shows  that  the  sediment  characteristics  do  change 
along  these  beaches,  and  differences  in  offshore  beach  slopes  between 
Wrightsville  Beach  and  Carolina  Beach  were  discussed  in  Section  II.  Th 
second  possible  conclusion,  and  probably  the  more  dominant  one  for  this 
study,  is  that  the  value  of  /9  is  very  sensitive  to  the  method  of  com¬ 
putation  of  the  variables  in  'the  rates  q^/CdP^/dX).  In  particular, 
errors  inherent  within  the  refraction  analysis  technique  can  result  in 
significant  spatial  variation  of  the  energy  flux  and  hence  in  the 
dP^/dX  values.  This  variation  is  then  reflected  in  the  spatial 
variation  of  the  /9  values. 


Table  16.  Values  of  /3  for  Wrightsville  and  Carolina  Beaches. 


Beach 

Values 

of  /3  in  units  of  m^-s/N-yr 

Best  fit 

Lower  bound 

Wrightsville 

300 

200 

400 

Carolina 

900 

300 

ANNUAL  VOLUME  CHANGE  PER  UNIT  LENGTH  OF  BEACH  <m3/yr/m) 


.  COMPUTED  VALUES  USING#  >1500 
.COMPUTED  VALUES  USING#  >900 
.  COMPUTED  VALUES  USING  #  -  300 
,  MEASURED  DATA  VALUES 


NOTE:  UNITS  OF 


ter] 


DISTANCE  ALONG  BEACH  (km) 

Figure  39.  Comparison  of  measured  and  computed  volumetric 
change  along  Wrightsville  Beach. 


.  COMPUTED  VALUES  USING  0 

—  — —  COMPUTED  VALUES  USING  /9 

—  —  —  —  —  COMPUTED  VALUES  USING  /3 

*  »  *  MEASURED  DATA  VALUES 


=  400 
=  300 
*=200 


Figure  40.  Comparison  of  measured  and  computed  volumetric 
change  along  Carolina  Beach. 
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Comparison  of  the  results  of  this  study  with  those  of  Jarrett's 
(1977)  are  encouraging.  Although  Jarret  calculated  his  /3  value  based 
only  on  the  midsection  of  Wrightsville  Beach,  his  value  of 
/3  •  418  m3-s/N-yr  is  approximately  equal  to  the  upper  limit  of  the 
value  of  /9  for  Wrightsville  Beach  as  predicted  by  this  study. 

5 .  Sediment  Budget. 

To  illustrate  the  application  of  the  sediment  transport  model  in 
estimating  the  northerly  and  southerly  longshore  transport  rates  and  the 
quantity  of  material  lost  into  the  adjacent  inlets,  sediment  budgets 
using  littoral  cells  of  finite  length  along  Wrightsville  and  Carolina 
Beaches  were  performed.  Each  beach  was  divided  into  the  three  cells 
which,  as  described  in  Section  IV,  4,  best  represent  the  long-term 
volumetric  changes  along  those  beaches.  Losses  from  the  active  profile 
due  to  a  rise  in  sea  level,  losses  from  the  beach  due  to  inlet  trapping, 
and  losses  or  gains  in  each  cell  due  to  longshore  sediment  transport 
were  all  considered.  The  long-term  excursion  rates  which  were  used  to 
determine  the  annual  volumetric  beach  change  for  each  cell  were 
calculated  by  eliminating  identified  excursions,  both  within  the  project 
boundaries  and  along  downdrift  beaches,  due  to  the  placement  and 
subsequent  initial  erosion  of  beach  fills.  Consequently,  the 
contributions  to,  and  the  commensurate  offshore  losses  from,  the  overall 
sediment  budget  due  to  beach-fill  operations  were  addressed  and  do  not 
need  to  be  further  incorporated  into  the  sediment  budget  equations. 
Aeolian  losses  were  considered  inconsequential  (U.S.  Army  Engineer 
District,  Wilmington,  1977)  and  also  were  not  included.  An  inherent 
assumption  within  this  approach  to  developing  a  sediment  budget  is  that 
offshore  losses  due  to  ongoing  sorting  of  freshly  exposed  beach  face  is 
minimal.  This  assumption  is  addressed  later  in  Section  VI  and  was  found 
to  be  valid. 

Based  on  the  concept  of  maintenance  of  an  equilibrium  profile  under 
rising  sea  conditions  (Bruun,  1962),  the  annual  volumetric  loss  of 
sediment  due  to  a  sea  level  rise  is  shown  in  Tables  6  and  17.  Losses 
due  to  wave  overtopping  occurred  only  along  the  northern  section  of 
Carolina  Beach.  Aerial  photos  taken  in  May  1964  and  November  1974  were 
used  to  estimate  the  bayward  excursion  of  the  bayside  shoreline. 

Results  from  that  analysis  indicated  that  approximately  4,600  ro3/yr 
was  lost  from  the  Oceanside  of  Carolina  Beach  (U.S.  Army  Engineer 
District,  Wilmington,  1977). 


Table  17.  Annual  volumetric  changes  in  beach-cell  volume  and  losses 


due  to  sea  level  rise  and  wave  overtopping. 


Beach 
cel  l 

Change  in  beach¬ 
cell  volume 
(m3/ yr ) 

Loss  due  to 
sea  level  rise 
(m3/ yr ) 

Loss  due  to 
wave  overtopping 
(m3/yr ) 

Wrightsville  (north) 

-24,430 

2,289 

_  __ 

Wrightsville  (central) 

-. 7,530 

1,873 

— 

Wrightsville  (south) 

-12,370 

1,457 

— 

Carolina  (north) 

♦104,500 

700 

4,600 

Carolina  (central) 

-269,750 

1,582 

— 

Carolina  (south) 

-107,970 

2,285 

— 

The  sediment  budget  equations  for  a  typical  beach  cell  (see  Fig.  41) 

are: 

Sediment  sources:  Qn-l,n  +  Qn+1 ,n 

Sediment  losses:  Qn,n-1  +  Qn,n+1  +  SLn  +  0Tn 

Annual  volumetric  beach  change: 

Vn  ~  Qn-l,n  +  Qn+l,n  “  Qn,n-1  "  Qn,n+1  _  JjLn  ~  ®Tn  (13) 

where  n,  n-1,  and  n+1  are  individual  beach  cells,  SLn  is  the  annual 
sediment  loss  from  cell  n  due  to  the  rise  in  sea  level,  0Tn  is  the 
annual  sediment  loss  from  cell  n  due  to  wave  overtopping,  and  Qn  n+j 
is  the  annual  longshore  sediment  transport  from  cell  n  into  cell  n+1. 

Equation  (11)  is  used  to  predict  the  quantity  Q  between  littoral 
cells  located  on  a  continuous  beach;  however,  a  problem  with  this 
formulation  arises  when  a  cell  boundary  borders  an  inlet,  weir  jetty, 
headland,  etc.  In  these  situations,  the  actual  quantity  of  sediment 
moving  in  the  littoral  drift  may  be  less  than  that  predicted  by 
equation  (11)  and  so  a  modification  must  be  incorporated  into  the 
sediment  budget  equations.  The  actual  longshore  sediment  transport 
rate,  Qa,  is  related  to  the  potential  longshore  sediment  transport  rate 
by  the  "efficiency  factor,"  a,  such  that 

Qa  =  a  (/SPx)  (14) 

Along  straight  and  continuous  beaches,  the  value  of  a  must  be  unity, 
however,  at  inlets  and  other  sediment  traps,  its  value  is  less  than  or 
equal  to  one.  In  extreme  cases  of  total  sediment  removal,  the  value  of  a 
is  zero.  The  solution  of  all  sediment  budget  equations  for  s  set  of 
littoral  cells  defines  the  values  of  a  at  each  cell  boundary. 

The  sediment  budget  schemat izat ions  for  Wrightsville  and  Carolina 
Beaches  are  shown  in  Figure  42.  The  values  of  the  northerly  and 
southerly  components  of  the  longshore  energy  flux  at  each  littoral  cell 
boundary  are  shown  in  Table  18.  The  values  of  /9  used  in  the  longshore 
sediment  transport  equations  were  j8  =300  for  Wrightsville  Beach  and 
/3  =900  for  Carolina  Beach.  The  measured  volumetric  change  within  each 
cell,  the  annual  volumetric  loss  due  to  sea  level  rise,  and  the  loss  due 
to  wave  overtopping  are  shown  in  Table  17. 

The  sets  of  a  values  at  each  inlet  boundary  (i.e.,  a  j  2  and 
a2  l»a4,5  an<l  a5,4;  and  07,8  andag(7)  cannot  be  uniquely 
determined  (there  are  more  unknowns  than  equations)  and  therefore,  the 
values  of  one  efficiency  factor  of  each  pair  must  be  assumed.  For  an 
unimproved  inlet  (i.e.,  no  jetties,  weirs,  etc.),  it  was  assumed  that 
all  sediment  contained  within  the  littoral  drift  system  entered  the 
inlet  cell.  In  this  case,  the  northerly  longshore  transport  from  the 
northern  ends  of  Wrightsville  and  Carolina  Beaches  was  assixned  to  enter 
Mason  and  Carolina  Beach  Inlets,  respectively.  Consequently,  a  2,1 
and  ag^  were  set  equal  to  one  and  the  sediment  budget  equations 
solved  resulting  in  the  values  ofa^2'*0*09  and  a  7  t8=0 . 31 . 
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Note:  Arrows  indicate  direction  of  sediment  movement. 


Figure  42.  Sediment  budget  schematics  for  Wrightsvllle  Beach 
and  Carolina  Beach. 


Table  18.  Energy  flux  values  at  cell  boindaries. 


Beach  cell 

Cell 

Cell 

ii'ffrrmoaiai: 

No. 

boundaries 

(Ian) 

Notation 

►tagnitule 

(N-m/s/m) 

Notation 

Mtgnitide 

(N-m/s/m) 

Northern  boundary 
(Mason  Inlet) 

1 

0.0 

— 

— 

Pl,2 

421.5 

Wrightsville  (north) 

H 

0.0-2. 5 

P2,l 

115.8 

P2,3 

287.6 

Wrightsville  (central) 

B 

2.5-4 .8 

p3,2 

290.0 

p3,4 

424.3 

Wrightsville  (south) 

B 

4.8-6. 7 

p4,3 

175.3 

p4,5 

446.3 

Masonboro  Inlet 

5 

6. 7-7.2 

P5,4 

352.8 

— 

— 

Carolina  Beach  Inlet 

7 

19.7-20.5 

— 

— 

^,8 

310.2 

Carolina  (north) 

8 

20.5-21.5 

p8,7 

226.7 

p8,9 

165.6 

Carolina  (central) 

21.5-24.3 

p9 ,8 

418.6 

p9,10 

350.0 

Carolina  (south) 

t! 

24.3-27.3 

P 

r  10 ,9 

305.0 

p 

10,11 

386.5 

Kure  and  Fort  Fisher 

n 

27.3-42.0 

e 

Pa 

224.2 

These  values  indicate  that  approximately  90  percent  of  the  potential 
southerly  longshore  sediment  transport  remained  trapped  in  Mason  Inlet 
and  70  percent  remained  in  Carolina  Beach  Inlet. 

The  north  jetty  at  Masonboro  Inlet  was  completed  in  spring  1966  and 
consists  of  a  rubble-mound  outer  section  and  a  low  concrete  sheet-pile 
inner  or  weir  section.  The  design  of  this  weir  jetty  and  the  dredging 
of  material  from  the  deposition  basin  on  the  inlet  side  of  the  weir  have 
caused  a  reduction  in  the  northward  sediment  bypassing  to  near  zero 
(U.S.  Army  Engineer  District,  Wilmington,  1977).  Therefore,  a  5  4 
was  set  equal  to  zero  and  the  solution  of  the  sediment  budget  equations 
gave  a 4  5=0.64.  This  means  that  approximately  two-thirds  of  the 
potential  littoral  drift  passes  over  or  around  the  weir  jetty  into 
Masonboro  Inlet  and  one-third  remains  trapped  on  the  southern  end  of 
Wrightsville  Beach,  providing  a  source  of  material  for  northerly 
transport.  Table  19  gives  the  a  values  for  the  Wrightsville  Beach  and 
Carolina  Beach  sediment  budgets. 


Table  19.  Efficiency  factors  a  for  Wrightsville  and  Carolina  Beach 
sediment  budgets. 


Beach  cell 

Cell 
No . 

[  Southerly  transport 

Notation 

Factor 

Factor 

Mason  Inlet 

1 

— 

— 

1,2 

Wright mril  le  (north) 

H 

2,1 

1.0 

2,3 

Wrightsville  (central) 

B 

3.2 

1.0 

3,4 

1.0 

Wrightsville  (south) 

B 

4,3 

1.0 

4,5 

0.64 

Masonboro  Inlet 

■> 

5,4 

1 

- 

Carolina  Beach  Inlet 

7 

- 

- 

1  7,8 

0.31 

Carolina  (north) 

8 

8,7 

1.0 

!  8.9 

1.0 

Carolina  (central) 

9 

1  9,8 

1.0 

i 

1  O 

Carolina  (south) 

to 

if.” 

1.0 

1  ’  j.  1 1 

;  0 

Kure  and  Fort  Fisher 

11 

I  u.n 

1.0 

— 

— 
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A,:.:  were  performed  to  include  Masonboro,  Kure,  and  Fort  Fisher 

Beaches  into  one  continuous  sediment  budget  analysis;  however,  the  lack 
of  reliable  long-term  volumetric  change  data  along  those  beaches  meant 
that  large  and  somewhat  arbitrary  changes  in  either  the  volumetric 
excursion  rates,  energy  flux  values,  or  /3  values  were  needed  to  balance 
all  sedinunt  budget  equations.  Because  of  these  changes,  the  results 
were  not  meaningful  and  are  not  presented. 


VI.  BEACH-FILL  PERFORMANCE 

All  beach  tills  placed  along  the  study  area  between  1965  and  1975 
were  discussed  in  Section  II;  Table  4  and  Figure  8  of  Section  IV  show 
additional  detailed  information  on  their  location  and  time  of  placement. 
The  beach  rills  are  also  discussed  in  Vallianos  (1970),  U.S.  Army 
Engineer  District,  Wilmington  (1970,  1974,  1977),  and  Jarrett  (1977). 
Information  presented  in  this  section  is  based  on  the  quantitative 
interpretation  of  the  excursion  distance  analyses  of  the  197C  beach  fill 
on  Wrightsville  Beach  and  of  the  1965  and  1971  beach  fills  on  Carolina 
Beach.  There  was  insufficient  repetitive  profile  information  for  the 
other  fills  to  allow  excursion  analysis  and  subsequent  fill  performance 
evaluation. 

The  1970  beach  fill  along  the  central  part  of  Wrightsville  Beach  was 
the  best  documented  (in  terms  of  repetitive  beach  surveys  before  and 
after  placement  of  fill  material)  beach-fill  project,  and  the  excursion 
distance  plots  of  profiles  WB13  to  WB29  (App.  A)  show  the  response  of 
the  beach  lo  this  fill.  Sequential  profiles  showing  the  post-fill 
behavior  at  profile  WB-15  are  presented  in  Appendix  F.  All  relevant 
data  from  all  of  these  plots  are  summarized  in  Figure  43  which  shows  the 
spatial  variation  along  the  beach  of  the  initial  fill  excursion,  the 
percent  total  initial  losses,  the  net  excursion  after  initial  losses, 
the  long-term  erosion  rate,  and  the  value  of  the  exponential  decay 
constant,  k.  All  values  in  the  figure  are  averaged  from  the  MLW,  MSL, 
and  MHW  excursion  distance  plot  of  each  profile  located  along  the 
central  section  of  Wrightsville  Beach. 

The  average  initial  fill  excursion,  as  defined  by  the  first  measure¬ 
ments  taken  efter  fill  placement,  was  76.6  meters,  and  the  distribution 
of  the  fill  along  the  beach  was  almost  triangular.  The  maximum  initial 
excursion  was  approximately  125  meters  in  the  middle  and  the  excursion 
at  the  project  boundaries  was  approximately  50  meters.  Figure  43  shows 
that  beach  excursions  were  measurable  along  the  beaches  on  either  side 
of  the  project  boundaries  soon  after  the  initial  fill  placement.  These 
edge  excursions  indicate  that  some  of  the  material  placed  within  the 
project  limits  of  the  fill  quickly  spread  laterally  to  the  adjacent 
beaches.  The  average  fill  excursion  remaining  on  the  beach  face,  after 
all  initial  losses  had  occurred  (approximately  2  years),  was  15.5  meters 
with  a  maximum  retention  of  29  meters  in  the  middle  of  the  fill.  This 
means  that  80  percent  of  the  initial  fill  was  lost  due  to  sorting,  slope 
readjustment,  and  lateral  spreading.  The  southern  end  of  the  fill 
exper  i  .-.i,  ed  tiie  highest  initial  loss  of  90  percent  where  only  5  meters 
"f  •  <  ii  1 1  remained  after  approximately  1.5  years. 
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INITIAL  FILL  EXCURSION 


k 


Wrightsville  Beach  1970  beach  fill 


' 


During  the  calculation  of  volumetric  change  between  two  subsequent 
profiles,  based  upon  the  application  of  the  volumetric  equivalent  factor 
to  the  MLW-MHW  contour  excursions  of  those  profiles,  an  assunption  of 
self  similarity  in  profile  shape  was  employed.  In  other  words, 
volumetric  changes  were  assumed  to  occur  only  as  a  result  of  horizontal 
displacement  of  the  profile  and  not  to  the  redistribution  of  material 
from  the  upper  beach  face  offshore,  a  phenomenon  which  occurs  during  the 
slope  readjustment  phase  of  the  beach-fill  response.  Consequently,  the 
total  initial  volumetric  loss  for  the  fill  may  be  slightly  less  than  the 
80  percent  value;  however,  the  average  initial  loss  in  beach  face 
position  is  still  80  percent  of  the  fill  excursion. 

The  adjustment  during  the  design  phase  of  the  project  for  the 
expected  sorting  losses  was  accomplished  by  applying  a  factor  known  as 
the  critical  ratio  (or  beach-fill  factor)  to  the  required  volume  of 
beach  fill.  The  critical  ratio  is  simply  an  estimate  of  the  quantity  of 
borrow  material  required  to  yield  1  cubic  meter  of  beach  material  having 
granulometric  characteristics  similar  to  the  native  beach.  The  value 
calculated  for  the  Banks  Channel  borrow  site,  and  which  was  applied  to 
the  Shell  Island  borrow  material,  was  2.5  (U.S.  Army  Engineer  District, 
Wilmington  1977).  This  means  that  2.5  cubic  meters  of  fill  material  was 
required  to  produce  1  cubic  meter  of  fill  material  on  the  beach  after 
sorting;  i.e.,  a  60-percent  sorting  loss  was  expected. 

A  modification  to  the  original  fill-factor  formulation  was  developed 
by  James  (1965)  and  has  now  been  incorporated  into  modern  beach-fill 
design  practices  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  1977).  Granulometric  data  from  profiles  taken  in  July 
1969  just  before  the  fill  and  samples  taken  from  profiles  along  the  fill 
just  after  placement  are  shown  in  Table  20.  These  values  were  used  to 
calculate  the  adjusted  fill  factor,  RA,  from  Figure  5.3  of  the  SPM. 

The  value  of  the  adjusted  fill  factor  was  RA=3.0,  which  implies  that 
66  percent  of  the  initial  fill  was  lost  to  sorting.  The  new  adjusted 
fill  factor  predicted  larger  sorting  losses  than  did  the  older  formu¬ 
lation;  however,  both  methods  predicted  losses  that  were  lower  than  that 
measured.  Assuming  that  these  formulations  are  correct,  then  losses  in 
addition  to  sorting  (slope  readjustment,  lateral  spreading,  etc.),  are 
significant  and  must  be  included  in  the  beach-fill  design. 


Table  20.  Granulometric  data  for  Wrightsville  Beach  1970  beach  fill. 


Granulometric 

conditions 

Date 

Profile 

Compos i te 
mean  grain 
s  ize  H 

(in  phi  units) 

Compos ite 

standard 
deviation  s 
(in  phi  units) 

Before  fill 

July  1969 

WB16 

1.53 

0.41 

WB29 

1.52 

0.83 

Prefill  coi 

uposite  values 

1.52 

0.87 

After  fill 

— 

Aug.  1970 

WB13 

2.23 

0.49 

WB19 

1.64 

0.76 

WB25 

1.78 

0.83 

Prefill  composite  values 

1.88 

0.69 

The  values  in  Table  20  were  also  used  to  calculate  James'  (1974) 
renourishment  factor,  Rj=1.9.  This  factor  expresses  the  ratio  of 
the  retreat  rate  of  the  beach  after  fill  placement  to  the  retreat  rate 
before  beach-fill  operations.  However,  in  its  derivation,  James  (1974) 
assumed  that  the  postfill  retreat  rate  was  linear  and  not  exponential. 
Therefore,  its  value  cannot  be  compared  to  the  results  of  this  study. 

The  relative  changes  in  the  upper  beach-face  angle  (from  MHW  to  MLW) 
after  fill  placement  were  measured  for  profile  WB17.  Figure  44  shows 
that  immediately  after  placement  the  average  beach  face  angle  was  1  on 
57,  which  was  flatter  than  the  prefill  angle  of  1  on  35.  The  beach 
angle  changed  fairly  rapidly  during  the  first  6  months  after  placement, 
and  after  9  to  12  months,  the  difference  in  the  average  beach  angle  at 
that  time  and  the  long-term  beach  face  angle  was  less  than  the  expected 
difference  due  to  seasonal  fluctuations.  It  is  apparent  that  a  signifi¬ 
cant  proportion  of  the  upper  beach  slope  adjustments  and  sorting  losses 
occurred  during  the  first  9  to  12  months.  After  that  period,  the  upper 
beach  face  retreated  with  a  fairly  constant  slope. 

The  value  of  the  exponential  decay  constant,  determined  from  the 
average  of  the  individual  k  values  for  each  of  the  MLW,  MSL,  and  MHW 
excursion  plots  from  each  profile,  was  k=0.66.  Substituting  this  value 
into  equation  (4),  together  with  4=0.8  and  St=0.954fi>  gave  t^=1.8  years 
i.e.,  effectively  all  initial  losses  due  to  sorting,  slope  adjustment, 
and  lateral  spreading  occurred  during  the  first  1.8  years  after  fill 
completion.  Substituting  k=0 . 66 ,  £=0 . 8 ,  and  Et=0  into  equation  (6) 
produces  t=4.06  years.  This  means  that  the  beach  face  eroded  back  to 
its  original  prefill  position  4.06  years  after  fill  completion,  and  that 
the  beach-fill  project  effectively  "bought"  this  time  for  the  beach 
segment  within  the  project  boundaries  by  artificially  placing  sand  on 
the  beach.  This  is  in  agreement  with  observed  behavior.  Between 
October  1970  and  December  1974,  an  estimated  91  percent  of  the  initial 
beach  fill  was  lost  (U.S.  Array  Engineer  District,  Wilmington,  1977),  and 
the  sequential  beach  profiles  in  Appendix  F  show  that  by  April  1974  the 
location  of  profile  WB15  was  approximately  in  its  pre-1970  beach-fill 
position.  Only  a  few  percent  of  the  initial  fill  was  retained  above  the 
MHW  contour  after  4  years  and,  unfortunately,  little  information  is 
available  to  describe  the  changes  in  offshore  bathymetry.  Downdrift 
beaches  benefited  from  the  fill  due  to  alongshore  transport  away  from 
the  fill  site.  However,  quantification  of  this  benefit  was  not  possible 
due  to  the  masking  effect  of  the  seasonal  variations  in  beach  position. 

Assuming  that  only  slope  and  sorting  adjustments  occurred  during  the 
first  9  to  12  months,  then  solving  equation  (4)  for  St  at  t  ^ =0 . 7  5 
and  t2=1.0  indicates  that  54  to  62  percent  of  the  total  initial  fill 
volune  was  lost  to  sorting  and  slope  adjustment.  This  range  compares 
favorably  with  the  values  of  60  to  66  percent  sorting  loss  estimated  by 
the  adjusted  fill  factor  and  critical  ratio  techniques,  respectively. 

The  rate  of  initial  loss  of  beach  material  was  not  constant  along 
the  length  of  the  beach-fill  project.  The  k  values  calculated  for 
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(1  On  37) 


Response  of  foreshore  slope  offer  Wrightsville  Beach  1970  beach 


profiles  near  the  ends  of  the  fill  tended  to  be  slightly  high.-i  *uuu 
those  for  profiles  located  in  the  middle.  This  implies  that  tie-  •  m  f 
the  fill  eroded  at  a  slightly  faster  rate  than  did  the  center  ,  wl  h  an 
be  expected  since  the  relative  changes  in  beach  angle  and  nearsn  i.- 
bathymetry  at  the  ends  are  greater  than  the  relative  changes  m  tie 
center,  and  thus  cause  greater  concentration  of  wave  energy  arid  sediment 
transport.  Together  with  the  fact  that  20-  to  30-meter  excursions 
occurred  on  either  side  of  the  fill  soon  after  placement,  this  inlotnia- 
tion  supports  the  concept  that  significant  quantities  of  fill  materia', 
spread  laterally  from  the  fill  ends.  It  should  be  acted,  however  ,  that 
nonhomogeneity  in  the  fill  material  properties  may  have  been  the  real 
cause  of  the  variation  in  the  rate  of  initial  loss  along  the  project 
length.  Approximately  70  percent  of  the  fill  material  was  obtained  from 
a  shoal  in  the  Banks  Channel,  and  the  balance  which  was  extremely  fine 
sand  of  poor  beach-fill  quality  was  obtained  from  the  sound  area  behind 
Shell  Island  (U.S.  Army  Engineer  District,  Wilmington,  1977). 

The  most  significant  feature  of  the  variation  in  long-term  excursion 
rate  along  central  Wrightsville  Beach  is  that  the  rate  calculated  for 
the  1963  to  1975  decade  (i.e.,  5  years  before  and  5  years  after  fill 
placement)  was  significantly  higher  in  the  vicinity  of  the  fill  than 
along  adjacent  beach  sections.  This  means  that  the  reason  for  the  high 
erosion  rates,  which  existed  before  and  probably  resulted  in  the  need  of 
the  1970  fill,  still  existed  after  1970  and  caused  high  annual  sediment 
losses  to  the  fill. 

There  are  two  possible  causes  for  these  localized  higher  erosion 
rates.  In  1965,  the  north  jetty  of  Masonboro  Inlet  was  completed  and 
effectively  cut  all  northward  sand  transfer  from  Masonboro  Island  to 
Wrightsville  Beach.  Consequently,  Wrightsville  Beach  suffered  higher 
erosional  losses  since  1965  due  to  the  partial  lack  of  sediment  supply. 
South  of  the  fill  the  growth  of  the  accretion  fillet  may  have  offset  t  lie 
increased  erosional  trends;  however,  the  same  is  not  true  for  the  area 
adjacent  to  the  north  fill  boundary. 

An  oblique  aerial  photo  of  Wrightsville  Beach  taken  between  1968  and 
1969  (Fig.  45),  shows  a  significant  deviation  in  the  present-day  shore¬ 
line  alinement  near  the  center  of  the  island.  The  uniform-width  dai k 
band  between  the  beach  and  the  seawardmost  houses  is  the  grassed  part  of 
the  constructed  dune  of  the  1965  beach-fill  project.  The  mi  sal inement 
of  the  north  end  of  the  Wrightsville  Beach  fill,  relative  to  the  present 
tendency  of  the  shoreline,  resulted  from  Moore  Inlet  which,  prior  t  >  its 
artificial  closure  in  1965  as  part  of  the  overall  beach  nour i shmen t 
plan,  was  located  just  north  of  arrow  A.  The  closure  of  Moore  Inlet 
eliminated  the  interaction  between  tidal  and  littoral  forces  in  this 
area,  which  had  existed  since  1887  and  which  had  combined  to  form  a 
seaward  concavity  in  the  shoreline  alinement  immediately  south  of  the 
inlet.  Erosion  prior  to  the  1965  beach  fill  exposed  the  northern 
building  line  of  the  township  of  Wrightsville  Beach  and  so  the  aline¬ 
ment  of  the  1965  beach  fill  was  forced  to  follow  this  line,  thus  causing 
a  bulge  in  the  resulting  beach  planform.  Arrow  B  points  to  profile 
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Figure  45.  View  of  Wrlghtsvllle  Beach  looking  north-northeast. 


WB36,  the  approximate  start  of  the  alinement  problem.  Between  1965  and 
1970,  the  beach  on  the  north  side  of  the  Masonboro  Inlet  jetty  accreted, 
however,  the  central  island  bulge  and  alinement  problem  remained.  The 
1970  beach  fill  was  placed  approximately  between  arrows  A  and  A'  ,  thus 
reinforcing  the  beach  alinement  problems.  The  greater  relative  change 
in  beach  planform  and  nearshore  bathymetry  in  the  central  section  of  the 
island  from  1965  to  1975  resulted  in  higher  wave  activity  and  erosionai 
trends . 

Natural  beach  processes  tend  to  focus  on  and  smooth  out  irregular¬ 
ities,  thus  creating  a  smoothly  curving  beach  as  is  idealized  by  the 
dashline  in  Figure  45.  The  high  rates  of  erosion  and  initial  losses 
associated  with  the  1970  beach  fill  may  not  be  typical  of  all  beach 
fills,  but  may  have  been  partly  caused  by  the  exposure  to  increased  wave 
attack  due  to  the  misalineraent  of  the  beach  planform.  The  resulting 
implication  means  that  if  improvement  in  performance  of  a  future 
beach  fill  located  in  the  same  area  is  desired,  then  additional  fill 
should  be  placed  along  the  adjacent  beaches,  as  shown  by  the  dot-dash 
line  in  Figure  45,  to  remove  the  alinement  problem.  This,  however,  may 
not  be  an  economically  feasible  solution. 

Information  obtained  from  the  postfill  beach  response  was  used  to 
examine  the  assumption  in  the  sediment  budget  analysis  that  offshore 
losses  due  to  sorting  of  freshly  exposed  beach  material  were  minor. 
Equation  (6)  showed  that  4.06  years  after  the  fill  placement,  the  beach 
returned  to  its  prefill  position  and  with  approximately  the  same  near¬ 
shore  profile  (Fig.  44).  This  means  that  whatever  came  into  the  fill 
area  during  the  4.06-year  period  was  transported  out  by  the  end  of  that 
time . 

The  sources  of  sediment  include  longshore  transport  into  the  fill 
region,  material  placed  during  the  beach-fill  operations,  and  material 
brought  ashore  by  seasonal  onshore  transport.  Losses  of  sediment 
include  longshore  transport  out  of  the  fill  region,  losses  due  to 
sorting  of  the  beach  fill,  seasonal  losses  due  to  offshore  transport, 
losses  due  to  the  rising  sea  level,  overwash,  and  aeolian  processes. 
Since  the  pre-  and  end-of-period  profiles  had  approximately  the  same 
shape,  the  net  volimetric  changes  due  to  slope  readjustment  were  zero. 
Over  an  even  4-year  period,  seasonal  changes  should  approximately 
balance  out,  and  so  within  the  limits  of  accuracy  of  this  study,  the  net 
on/offshore  contribution  was  set  to  zero.  Volumetric  gains  from  the 
beach  fill  (BF)  were  determined  from  surveys,  and  associated  sorting 
losses  (sorting)  were  calculated  using  the  adjusted  fill  factor 
(R^).  Losses  due  to  sea  level  (SL)  were  calculated  by  use  of 
Bruun's  (1962)  formulation.  Aeolian  and  washover  losses  were  near  zero. 
Since  the  net  volume  change  at  the  end  of  the  4.06-year  period  was  zero, 
then  the  net  volumetric  change  due  to  alongshore  transport  of  the 
boundaries  (Qin~Qout^  ®ust  equal  the  difference  between  these 
identified  sources  and  sinks  since  the  fill  area  was  away  from  active 
inlets,  jetties,  etc.;  i.e., 

Qin  “  Qout  +  BF  "  sorting  -  SL  *  0 
Qin  “  Qout  +  bf/rA  -  SL  =  0 
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Substituting  in  values  from  Tables  4  and  6  and  with  R^=3,  the  annual 
volumetric  change  due  to  alongshore  transport  was 


From  equation  (11) 


and  from  Figure  38 


and  hence 


Qin  ~  Qout  =  ~86»125  m3/yr 


Qin  Qout  =  0  (Plin  Plout) 

P ] •  -Pi  =  -300  N-M/S-M 
1in  Aout 

/3  =  28  7  m3-s/Nyr 


Within  the  limits  of  accuracy  of  both  the  data  and  technical  analyses, 
this  simple  postfill  sediment  budget  determination,  where  all 
contributions  to  the  sediment  budget  were  quantified,  produced  a  value 
of  /9  which  was  very  close  to  that  calculated  earlier  (mean  value  of 
/3  =  300)  when  using  long-term  beach  response  characteristics  and  where 
the  losses  due  to  sorting  of  "freshly"  exposed  native  beach  material  by 
ongoing  erosion  was  assumed  to  be  small.  Since  the  calculated  values  of 
/3  are  similar  and  they  come  from  analyses  of  two  distinct  phases  of 
beach  response,  these  results  support  the  contention  that  ongoing 
sorting  losses  during  the  long-term  response  phase  are  minimal. 

Analysis  of  the  spatial  variation  of  the  beach  response  to  the  1965 
and  1971  beach  fills  along  Carolina  Beach  was  not  possible  because  of 
insufficient  profile  information.  Results  for  the  1965  fill,  as  shown 
by  the  beach  photos  in  Figure  46,  were  determined  from  MLW ,  MSL,  and  MHW 
excursion  distance  plots  for  profiles  CB106  and  CB107  which  were  less 
than  0.5  kilometer  apart.  Consecutive  profiles  at  CB97  were  used  to 
determine  the  response  to  the  1971  beach  fill.  The  average  exponentia 
decay  constant,  the  average  initial  fill  excursion,  and  the  average 
long-term  erosional  rate  are  given  in  Table  12.  Substituting  these  val¬ 
ues  into  equation  (4)  indicates  that  most  initial  losses  occurred  during 
the  first  1.5  to  2  years  following  both  fills,  in  agreement  with 
observed  behavior  (U.S.  Army  Engineer  District,  Wilmington,  1970). 

Using  the  values  contained  in  Table  12  and  assuming  £ =0 .8 ,  equation  (6) 
predicts  that  2.4  years  and  2.25  years  after  the  1965  and  1971  fill 
projects,  respectively,  the  beach  face  eroded  to  approximately  its 
original  prefill  position.  These  values  are  in  reasonable  agrennent 
with  recorded  observations  on  the  loss  of  beach  fill  during  the  2  years 
following  each  fill  (U.S.  Army  Engineer  District,  Wilmington,  1977). 
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data  taken  immediately  after  fill  placement  in  1965, 
years  later,  are  shown  in  Table  21.  These  data  were 
a  critical  ratio  of  2.1  for  the  fill  material,  and 
55  percent  volumetric  loss  due  to  sorting  (U.S.  Army 
,  Wilmington,  1970).  Results  from  profile  '  B106  tend 
ercent  of  the  initial  excursion  was  lost  during  the 
ars,  close  to  the  design  value.  The  adjusted  fill 
(1975)  renourishment  factor  were  evaluated  from  the 
e  found  to  be  R^=I .02  and  Rj=0.25,  respectively.  For 
beach-fill  data,  the  adjusted  fill-factor  techniques 
of  expected  sorting  loss  signi f icant 1 y  lower  than  both 


VII.  SUMMARY  AND  CONCLUSIONS 


During  the  period  from  1964  to  1975,  2,952  repetitive  beach  profiles 
were  recorded  at  241  stations  between  Wrightsville  Beach  and  Fort 
Fisher  Beach.  The  total  length  of  Wrightsville  and  Carolina  Beaches 
represented  only  32  percent  of  the  total  length  of  the  study  area,  but 
nearly  70  percent  of  all  beach  profile  stations  and  89  percent  of  the 
total  number  of  recorded  profiles  were  located  along  these  two  beaches. 
Of  the  nearly  3,000  profiles,  only  4  percent  extended  beyond  the  MLW 
position  to  approximatley  the  -10  meter  contour.  As  a  consequence, 
volumetric  changes  representative  of  actual  changes  occurring  between 
successive  surveys  could  not  be  calculated  by  simply  measuring  the 
change  in  area  under  the  measured  profile  curves  because  significant 
changes  occur  below  the  low  water  line. 

The  positions  of  the  MHW,  MSL,  MLW,  -1.83  meters  (-6  feet), 

-3.66  meters  (-12  feet),  and  -5.49  meters  (-18  feet)  contours  were 
plotted  relative  to  a  fixed  base  line,  for  all  profiles.  The  excursion 
distance  of  each  contour  between  successive  profiles  is  indicative  of 
volumetric  change,  the  magnitude  of  which  is  found  by  applying  a 
volumetric  equivalent  factor,  calculated  from  changes  in  area  under  some 
profiles  which  repetitively  extended  out  into  deeper  water,  to  the  mean 
excursion  distance  value.  Due  to  the  poor  spatial  and  temporal 
distribution  of  profiles  along  Masonboro ,  Kure,  and  Fort  Fisher  Beaches, 
only  profiles  from  Wrightsville  and  Carolina  Beaches  were  used  in  the 
analysis  of  beach  response  and  volumetric  changes  associated  with  storms 
and  manmade  influences.  The  results  indicate  that  the  average  seasonal 
changes  along  Wrightsville  and  Carolina  Beaches,  measured  24  and 
17  meters,  respectively,  were  significantly  larger  than  the  long-term 
loss  (erosion)  rate  for  1  year.  In  addition,  the  response  of  these 
beaches  to  storm-induced  erosion  or  beach-fill  placement  was,  in  many 
instances,  very  short  in  duration  and  therefore  difficult  to  identify  in 
many  of  the  excursion  plots  which  had  poor  temporal  resolution. 

Most  of  the  beach  profile  data  are  not  a  result  of  one  coordinated 
and  well-planned  study,  but  rather  from  several  independent  and  over¬ 
lapping  studies.  The  following  recommendations  on  the  distribution  of 
beach  profile  surveys  are  based  on  comparison  of  adjacent  profiles  and 
are  made  so  that  the  most  efficient  use  of  manpower  and  money  can  be 
incorporated  into  future  beach  studies. 

The  spatial  separation  of  profiles  should  be  in  the  range  of  0.5  to 
1.0  kilometers,  if  possible,  along  straight  or  smoothly  varying 
stretches  of  beach.  Profiles  should  be  spaced  closer  in  areas  of  abrupt 
changes  in  beach  planform  (e.g.,  inlets,  headlands,  etc.)  or  in  areas 
where  historic  observations  indicate  large  relative  changes  in  beach 
pos i t ion . 

The  profiles  must  be  measured  with  sufficient  frequency  so  that 
seasonal  fluctuations  and  longer  term  trends  can  be  identified  and 
separated.  To  accomplish  this,  some  stations  (e.g.,  every  fourth)  must 
be  surveyed  frequently,  no  more  than  1  or  2  months  apart,  and  the  inter¬ 
mediate  stations  should  be  profiled  at  least  twice  a  year  (surveyed  at 
the  same  times  each  year) . 
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Some  of  the  profiles  which  are  surveyed  frequently  must  be  surveyed 
out  beyond  the  MLW  position  to  approximately  the  position  of  the 
-10  meter  contour.  These  long  profiles  are  necessary  to  establish  the 
actual  volumetric  changes  for  the  entire  active  profile,  and  hence,  used 
to  calculate  the  volumetric  equivalent  factor  applied  to  the  intermediate 
profiles . 

If  the  seasonal  variation  in  beach  excursion  is  larger  than  the 
long-term  trends,  then  profile  data  must  be  collected  for  a  minimum  of  2 
to  3  years  for  both  processes  to  be  quantified.  Greater  variability  in 
the  data  necessitates  longer  collection  periods. 

For  projects  with  tight  budget  constraints,  a  few  profiles  located 
in  key  positions  and  surveyed  frequently  will  provide  a  better  data  base 
than  more  profiles  surveyed  infrequently. 

Wave  gage  data  collected  at  Johnnie  Mercer's  Pier  and  LEO  data  from 
Wrightsville  Beach  were  combined  to  develop  a  wave  climate  representa¬ 
tive  of  the  wave  conditions  found  along  the  study  area.  This  data  was 
refracted  in  to  shore  and  the  breaking  wave  conditions  were  used  to 
calculate  both  the  northerly  and  southerly  components  of  longshore 
energy  flux.  The  spatial  gradient  of  these  values  along  Wrightsville 
and  Carolina  Beaches  were  compared  with  the  long-term  (nonseasonal) 
volumetric  changes,  and  the  empirical  factor,  /3 ,  which  relates  the 
longshore  sediment  transport  rate  to  the  longshore  component  of  energy 
flux  was  calculated.  By  choosing  a  best-fit  value  of  )3  =300  and 
^3=900  m^-s/N-yr  for  Wrightsville  and  Carolina  Beaches,  respectively, 
plots  of  predicted  and  measured  volumetric  change  due  to  longshore 
sediment  transport  along  each  beach  showed  similar  trends,  although  the 
absolute  magnitude  at  any  beach  location  was  different. 

To  improve  the  accuracy  of  the  energy  flux  computation  in  future 
studies,  the  following  recommendations  on  desirable  refraction  model 
characteristics  should  be  utilized  or  developed. 

(a)  Variable  grid  cell  spacing  should  be  used  to  allow  coarse-sized 
computational  cells  in  deep  water  and  finer  cells  in  the 
nearshore  region  where  greater  relative  changes  in  bathymetry 
can  cause  instability  problems. 

(b)  The  effects  of  diffraction  and  tidal  currents  on  wave 
propagation  should  be  included. 

(c)  The  dynamic  interrelationship  between  both  the  nearshore 
bathymetry  and  shoreline  planform,  and  the  sediment  transport 
potential  of  the  incoming  waves  should  be  incorporated.  The 
present  static  boundary  condition  representat ion  of  the  shore¬ 
line,  used  in  refraction  analysis  programs,  does  not  allow  for 
any  change  in  shape  in  the  shoreline  due  to  increased  sediment 
transport  capabilities  as  a  result  of  increased  (focused)  wave 
activity.  Thus  changes  in  refraction  patterns  and  beach 
approach  angles  due  to  beach  response  between  different  sets  of 
wave  types  used  to  represent  seasonal  or  annual  conditions 
should  be  included. 
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Until  these  improvements  can  be  incorporated,  the  results  of  this 
study  indicate  that  the  additional  expenses  incurred  due  to  the  use  of  a 
large  number  of  wave  rays  and  high  resolution  in  the  bathymetric  data 
cannot  be  justified. 

Sediment  budgets  were  developed  for  Wrightsville  and  Carolina 
Beaches.  These  two  beaches  were  each  divided  into  three  littoral  cells 
in  which  the  response  of  the  beach  to  all  natural  and  man-influenced 
changes  was  fairly  similar.  Long-term  volumetric  changes  were  assumed 
to  be  the  result  of  differences  in  longshore  sediment  transport  rates, 
sediment  loss  to  wave  overtopping,  and  to  sea  level  rise.  Losses  due  to 
ongoing  sorting  of  beach  sediment  were  considered  minor.  Values  of  wave 
energy  flux  at  each  cell  boundary  were  multiplied  by  the  empirical 
factor  li  which  relates  the  longshore  transport  potential  to  the  long¬ 
shore  component  of  wave  energy  flux.  An  additional  efficiency  factor,  a 
which  relates  the  actual  volume  of  sediment  transported  to  the  potential 
amount  as  predicted  frou  the  energy  flux  analysis,  was  included  in  the 
sediment  budget  equations.  The  value  of  a  along  a  smooth  and  uninter¬ 
rupted  coastline  was  assumed  to  be  one  and  at  positions  where  a  coastal 
structure  (e.g.,  the  north  jetty  weir  at  Masonboro  Inlet)  or  where 
geologic  control  (availability  of  sediment  supply)  prohibit  transport, 
the  value  of  a  was  assumed  to  be  zero.  The  solution  of  the  sediment 
transport  equations  resulted  in  a  values  which  indicated  that  only 
two-thirds  of  the  gross  southerly  transport  along  Wrightsville  Beach 
spills  over  the  north  jetty  weir  into  Masonboro  Inlet  and  one-third  is 
either  trapped  along  the  southern  end  of  Wrightsville  Beach  or  locally 
transported  northward  by  wave  energy  reversals.  At  the  northern  ends  of 
Wrightsville  and  Carolina  Beaches,  only  10  and  31  percent  of  the 
potential  volume  of  sediment  is  transported  out  of  Mason  and  Carolina 
Beach  Inlets,  respectively.  If  better  volumetric  change  data  had  been 
available  for  Masonboro  Beach,  then  the  influence  of  Masonboro  and 
Carolina  Beach  Inlets  in  terms  of  their  inlet  trapping  potential  on  the 
supply  and  storage  of  sand  have  been  determined. 

Analyses  of  the  beach  prociles  taken  along  Wrightsville  Beach  after 
the  1970  beach  fill  indicate  several  components  of  beach  response.  The 
first  component  was  a  long-term  loss  rate  of  -3.8  meters  per  year  which 
was  approximately  equal  to  the  long-term  loss  rate  during  the  5-year 
period  prior  to  the  1970  fill  operations.  This  rate  was  much  higher  in 
the  immediate  vicinity  of  the  fill  than  along  adjacent  beaches  both 
during  the  prefill  and  postfill  periods,  and  indicated  that  the  fill 
placement  did  not  reduce  or  eliminate  the  problem  which  resulted  in  the 
need  for  a  fill,  but  rather  provided  recreational  opportunity  and  "bought- 
time"  for  the  properties  behind  the  project  boundaries. 

In  addition  to  the  long-term  component,  an  exponential  loss  of 
beach-fill  volume  was  recorded  during  the  first  1.5  to  2  years.  Excur¬ 
sion  plot  analysis  showed  that  about  80  percent  of  the  total  initial 
fill  was  eroded  during  this  period  of  rapid  initial  loss,  and  that 
severe  storm  erosion  was  not  the  primary  cause  for  the  very  high  initial 
loss  rate. 

The  first  set  of  profile  measurements  taken  after  fill  completion 
indicated  that  the  fill  material  was  placed  at  a  beach  angle  shallower 
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than  the  existing  1970  prefill  beach  slope.  During  the  first  9  to 
12  months  the  MHW-MLW  beach  slope  steepened  (and  retreated)  in  response 
to  the  seaward  sorting  of  fine  sand  grains  and  to  the  readjustment  of 
the  profile  slope  to  the  prevailing  wave  conditions.  After  this  period, 
the  upper  beach  face  retreated  with  only  minor  changes  in  beach  slope 
due  to  seasonal  wave  climate  influences. 

Sediment  characteristics  of  the  fill  and  native  beach  material  were 
used  to  calculate  a  value  for  the  adjusted  fill  factor  of  R^=3.0. 

This  value  indicates  that  66  percent  of  the  fill  material  can  be 
expected  to  be  lost  due  to  sorting;  however,  comparison  with  measured 
results  indicates  that  this  calculation  underestimates  the  initial  loss 
percentage.  In  addition  to  the  sorting  and  slope  readjustment  losses, 
significant  quantities  of  fill  material  were  lost  due  to  the  lateral 
spreading  of  material  onto  adjacent  beaches. 

An  oblique  aerial  photo  taken  before  the  1970  beach  fill  showed  that 
the  placement  of  the  fill  could  only  have  reinforced  the  beach  alinement 
problem  along  Wrightsville  Beach.  Since  1965,  the  beach  section  which 
suffered  the  localized  and  high  erosion  rates  protruded  from  the 
generally  smooth,  curving  beach  planform.  It  was  concluded  that  the 
relative  change  in  beach  planform  and  nearshore  bathymetry  resulted  in 
an  increase  in  localized  wave  activity,  sediment  transport  potential, 
and  erosional  trends,  and  that  this  phenomena  would  continue  until  the 
relative  change  in  beach  shape  is  eliminated.  Therefore,  it  appears 
that  the  continual  renourishment  of  this  section  perpetuated  the  problem 
of  increased  localized  wave  activity. 

This  study  showed  that  beach  losses  in  addition  to  the  expected 
losses  due  to  sorting  and  slope  readjustment  occurred  during  the  initial 
1.5-  to  2-year  response  phase.  It  appeared  that  lateral  spreading  of 
the  fill  material  onto  adjacent  beaches,  due  to  the  forced  protrusion  of 
the  beach  fill  out  beyond  the  general  beach  alinement,  resulted  in  these 
additional  significant  losses. 
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APPENDIX  G 


WAVE  REFRACTION  PLOTS 
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